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1. Introduction 



In the standard cosmological model, at very early times the Universe undergoes a quasi 
de Sitter exponential expansion driven by a scalar field, the infiaton, with an almost 
fiat potential. The quantum fiuctuations of this field are thought to be at the origin of 
both the Large Scale Structures and the Cosmic Microwave Background (CMB) fiuc- 
tuations that we are able to observe at the present epoch [1]. CMB measurements 
indicate that the primordial density fiuctuations are of order 10~^, have an almost 
scale-invariant power spectrum and are fairly consistent with Gaussianity and statistical 
isotropy [2, 3, 4, 5, 6, 7, 8, 9]. All of these features find a convincing explanation within 
the infiationary paradigm. Nevertheless, deviations from the basic single-(scalar)field 
slow-roll model of infiation are allowed by experimental data. On one hand, it is then 
important to search for observational signatures that can help discriminate among all 
the possible scenarios; on the other hand, it is important to understand what the theo- 
retical predictions are in this respect for the different models. 

Non-Gaussianity and statistical anisotropy are two powerful signatures. A random 
field is defined "Gaussian" if it is entirely described by its two-point function, higher 
order connected correlators being equal to zero. Primordial non- Gaussianity [10, 11] 
is theoretically predicted by infiation: it arises from the interactions of the infiaton 
with gravity and from self-interactions. However, it is observably too small in the 
single-field slow-roll scenario [12, 13, 152]. Alternatives to the latter have been 
proposed that predict higher levels of non-Gaussianity such as multifield scenarios 
[15, 16, 17, 18, 19, 20, 21], curvaton models [22, 23, 24, 25, 26, 27] and models with 
non-canonical Lagrangians [28, 29, 30, 31, 32]. Many efforts have been directed to the 
study of higher order (three and four-point) cosmological correlators in these models 
[33, 34, 35, 36, 37, 31, 152, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48] and towards 
improving the prediction for the two-point function, through quantum loop calculations 
[49, 50, 51, 13, 52, 53, 54, 55, 56]. From WMAP, the bounds on the bispectrum amplitude 
are given by -4 < < 80 [8] and by -125 < f^f < 435 [9] at 95% CL, respectively 
in the local and in the equilateral configurations. For the trispectrum, WMAP provides 
—5.6 X 10^ < qnl < 6.4 X 10^ [57] {qnl is the "local" trispectrum amplitude from cubic 
contributions), whereas from Large-Scale-Structures data —3.5 x 10^ < giyL < 8.2 x 10^ 
[58], at 95% CL. Planck [59] is expected to set further bounds on primordial non- 
Gaussianity. 

Statistical isotropy has always been considered one of the key features of the CMB 
fiuctuations. The appearance of some "anomalies" [60, 61, 62] in the observations 
though, after numerous and careful data analysis, suggests a possible a breaking of 
this symmetry that might have occurred at some point of the Universe history, possibly 
at very early times. This encouraged a series of attempts to model this event, preferably 
by incorporating it in theories of infiation. Let us shortly describe the above mentioned 
"anomalies". First of all, the large scale CMB quadrupole appears to be "too low" and 
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the octupole "too planar"; in addition to that, there seems to exist a preferred direction 
along which quadrupole and octupole are ahgned [63, 64, 60, 65, 66]. Also, a "cold spot" , 
i.e. a region of suppressed power, has been observed in the southern Galactic sky [61, 67]. 
Finally, an indication of asymmetry in the large-scale power spectrum and in higher- 
order correlation functions between the northern and the southern echptic hemispheres 
was found [68, 62, 69]. Possible explanations for these anomahes have been suggested 
such as improper foreground subtraction, WMAP systematics, statistical flukes; the 
possibilities of topological or cosmological origins for them have been proposed as well. 
Moreover, considering a power spectrum anisotropy due to the existence of a preferred 
spatial direction h and parametrized by a function g{k) as 



the five-year WMAP temperature data have been analyzed in order to find out what 
the magnitude and orientation of such an anisotropy could be. The magnitude has been 
found to be = 0.29 ± 0.031 and the orientation ahgned nearly along the echptic poles 
[70]. Similar results have been found in [71], where it is pointed out that the origin of 
such a signal is compatible with beam asymmetries (uncorrected in the maps) which 
should therefore be investigated before we can find out what the actual hmits on the 
primordial g are. 

Several fairly recent works have taken the direction of analysing the consequences, in 
terms of dynamics of the Universe and of cosmological fluctuations, of an anisotropic 
pre-inflationary or inflationary era. A cosmic no-hair conjecture exists according to 
which the presence of a cosmological constant at early times is expected to dilute any 
form of initial anisotropy [72]. This conjecture has been proven to be true for many (all 
Bianchi type cosmologies except for the the Bianchi type-lX, for which some restric- 
tions are needed to ensure the applicability of the theorem), but not all kinds of metrics 
and counterexamples exist in the literature [73, 74, 75]. Moreover, even in the event 
isotropization should occur, there is a chance that signatures from anisotropic inflation 
or from an anisotropic pre-inflationary era might still be visible today [76, 77, 78, 79]. 
In the same context of searching for models of the early Universe that might produce 
some anisotropy signatures at late time, new theories have been proposed such as spinor 
models [80, 81, 82, 83], higher p-forms [84, 85, 86, 87, 88, 89] and primordial vector field 



Within vector field models, higher order correlators had been computed in [90, 91, 92, 
93, 94] and, more recently, in [95, 96] for f/(l) vector fields. We considered SU{2) vector 
field models in [97, 98]. Non-Abelian theories offer a richer amount of predictions com- 
pared to the Abelian case. Indeed, self interactions provide extra contributions to the 
bispectrum and trispectrum of curvature fluctuations that are naturally absent in the 
Abelian case. We verified that these extra contributions can be equally important in a 
large subset of the parameter space of the theory and, in some case, can even become 
the dominant ones. 




(1.1) 



models. 
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The promising perspective of achieving more and more precise measurements for the 
cosmological observables thanks to Planck and future experiments and the search for 
signatures that may help identify the correct inflationary model, have also motivated 
studies of higher order corrections to cosmological correlation functions and to the power 
spectrum in particular. Indeed, loop corrections to the correlators arise from the interac- 
tions involving the fields during inflation and therefore carry some important information 
about the physics of the very early Universe. 

Loop corrections may lead to interesting effects which scale like the power of the 
number of e-folds between horizon exit of a given mode k and the end of infla- 
tion [99, 100, 101, 102, 103]. The interest in loop corrections to the correlators of 
cosmological perturbations generated during an early epoch of inflation has been re- 
cently stimulated by two papers of Weinberg [104, 105]. The reason is that one-loop 
corrections to the power spectrum of the curvature perturbation ( seem to show some 
infra-red divergences which scale like \n{kL), where is some infra-red comoving mo- 
mentum cut-off [106, 107, 108, 109]. However, it has been discussed in [110, 111] (see 
also [112]) that such potentially large corrections do not appear in quantities that are 
directly observable. 

As to the power spectrum of curvature perturbations, one-loop corrections have been 
computed in single-field slow-roll inflation by D. Seery [54, 55] and by N. Bartolo and 
myself [56], in single- field slow-roll inflation. In [56] we completed the analysis carried 
out in [54, 55], where the metric tensor fluctuations had been neglected for simplicity, by 
including them in the calculations and proving that their contribution is as important 
as the one from the scalar perturbations. In the context of loop-calculations, we have 
also been working on corrections to the power spectrum in theories with non-canonical 
Lagrangians, which allow for higher and possibly observable corrections [113]. 
It can be safely stated that in standard single-field slow-roll inflation, the perturbative 
expansion is well-behaved, in the sense that the agreement with observations found at 
tree-level for the power-spectrum is not spoiled by the radiative corrections and, on a 
more general basis, higher order loop corrections introduce smaller and smaller correc- 
tions as the perturbation series expansion progresses. This is not generically true in 
more general theories, such as for instance models with non-canonical Lagrangians, for 
which bounds need to be requested on the parameters of the theory in order to preserve 
the vahdity of the perturbative approach [114, 115]. 

This thesis collects the main results of our work on loop corrections to the power 
spectrum in theories of scalar inflation (Sees. 2 to 6) [56, 113], on anisotropic pre- 
inflationary cosmologies (Sec. 7) [78] and on primordial non-Gaussianity and anisotropy 
predictions from theories of inflation where vector fields can play a role in the production 
of the late time cosmological fluctuations (Sees. 8 to 12) [97, 98]. The (5N and the 
Schwinger-Keldysh formalisms are some of the main tools of our computation and will 
be briefly reviewed. 
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2. Schwinger-Keldysh formalism 



The temperature fluctuations in the CMB are rather small, of order 10~^. Theoretical 
predictions for the power-spectrum of curvature perturbations during inflation provide a 
very good match at tree level: this suggests that it is correct to use perturbation theory 
to evaluate cosmological correlators. A formalism conveniently employed to implement 
the perturbative approach is the Schwinger-Keldysh, also dubbed as "in-in" , formalism. 
It was first formulated in [49, 50, 51], later applied by J. Maldacena in [13] to the 
calculation of the bispectrum of curvature fluctuations and revived by S. Weinberg in 
[52, 53]. In this formalism the expectation value of a field operator 6(t) is given by 



{n\Q{t)\Q) 







f e 



,iJlHj{t')dt' 



0/W 



T e 



-iJlHi{t')dt' 







(2.1) 



where 1^2) represents the vacuum of the interacting theory, T and T are time-ordering 
and anti-ordering operators, the subscript / indicates the fields in interaction picture 
and Hj is the interaction Hamiltonian. The interaction picture has the advantage of 
allowing to deal with free fields only; the fields can be thus Fourier expanded in terms 
of quantum creation and annihilation operators 

6(f){x,t) = J Ske'''^ [a^Hkit) + a+g50fc(t) 

with commutation rules 

ag,a+] =(27r)V3)(^-fc"'0. 

The in-in formula has many similarities with the S-matrix in quantum field theory 
in terms of mathematical structure and perturbative approach, but they also have 
fundamental differences: the S-matrix corresponds to a transition amplitude between 
an initial and a final state; a cosmological correlation function is instead the expectation 
value of a given observable at a given time; moreover, asymptotic states in cosmology 
are only defined at very early times, when the same initial conditions as in Minkowsky 
spacetime apply for the free fields. 

Using the positive and negative path technique of the in-in formalism [52, 53], the 
expectation value above can be recast in the form 



(i]|e(t)|fi) 







T 







(2.2) 



where the plus and minus signs indicate modified Feynman propagators, i.e. modified 
rules of contraction between interacting fields; schematically we have 

(r(0i02---0n)) = [0i0j,0/0,„, ...], (2.3) 

ij,lm,... 

where the sum is taken over all of the possible sets of field contractions and 
0+(?7')0+(?7") = G> (77', 77") 6(77' - 77") + G<(7/',7/")e(7/" - r/'), 
G<(r/',r/"), 
G>(V,V'), 

G^{ri ,r] )Q{r] - r] ) + G^{r],7] )Q{'q - ?7 )• 



(p+{r]')(p-{r]"] 
(f)-{r]')(f)+{ri"] 



7 



In momentum space we have 

G>(V,V') = 50,(V)50KV'), 
G<{rl,ri')^5<i>l{rl)5<P,{rl'). 

It is important to remember that, when we apply this formahsm, the external fields are 
always supposed to be treated like +fields. 
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3. Scalar loop corrections to P<^ 



The power spectrum for the comoving curvature perturbation C is defined by 

{CkM,-,it)) = i2n)'P^{k)6^'\h + h) , (3.1) 

This and all other correlation functions presented in this thesis are computed using 
the SN formula, ({x) at a given time t can be interpreted as a geometrical quantity 
indicating the fiuctuations in the local expansion of the universe; in fact, if N{x,t*,t) 
is the number of e-foldings of expansion evaluated between times t* and t, where the 
initial hypersurface is chosen to be fiat and the final one is uniform density, we have 

C(f, t) = N{x, t*,t)- N{t*,t) = 6N{x, t). (3.2) 

The number of e-foldings N{x,t*,t) depends on all the fields and their perturbations 
on the initial slice. In principle, since the fields are governed by second order differen- 
tial equations, it should also depend on their first time derivatives, but if we assume 
that slow-roll conditions apply, then the time derivatives will not count as independent 
quantities. 

Let us apply Eq. (3.2) to the computation of in single-field slow-roll infiation (the 
Lagrangian for the scalar field is given by = {l/2)g^'^d^(f)du(f) — V{(l))) 



((E^..(.,nr),(i:^«..nr)). (3.3, 

The sums can be expanded to the desired order. Up to one loop we have 

(C.-;WC.-;(t))=iV^^^'(%-l%-^,)* 

+ i.iv(i)iv(3) J d'qdMScp,-^S<p,^cPpM,^_^^^.k,)* + ih ^ h) 

where a star indicates evaluation around the time of horizon crossing. Eq. (3.4) can 
finally be rewritten as [116, 109] 

(CkMd^)) = C^^fS^'Kh + h) [ (iV(l))' (Ptroc(A;i) + Ponc-loop(A:i) 
+ \ / 7?T3^t..ce(g)Pu.ee(|^l " <f\) 



2 V J J (2n) 

d^q 



+ N^"^ N^'^ P,Uk) j ^Ptree(g)] , (3.5) 
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Figure 1. Diagrammatic representation of the one loop corrections to the power 
spectrum of Scj} from scalar modes to leading (~ e°) order in slow-roll. 




Figure 2. Next-to-leading ^/e) order one loop corrections from scalar modes to 
the power spectrum of 6(j). 



PtTcc{k) is the tree level power spectrum (3.6) 



(3.6) 



where if* is the Hubble parameter evaluated at horizon exit (when k = aH). The 
variance per logarithmic interval in k is given by V{k) = [k^ /2'K'^)P{k). The one loop 
contribution to the power spectrum is given by 



where the first term on the right-hand side, Pscaiar, accounts for the contributions coming 
from the inflaton self-interactions and were computed by D. Seery in [108, 109] 



where gi and g2 are numerical factors. Their diagrammatic representation is given in 
Fig.l for the leading order and in Fig. 2 for the next-to-leading order corrections. The 
loop corrections Ptcnsor, arising from interactions between the tensor (graviton) modes 
and the scalar field, were ignored for simplicity in [108, 106], however they should be in- 
cluded since they are not slow-roll suppressed compared to loops of scalar modes. Their 
computation was presented for the first time in our paper [56] and will be reviewed in 
Sees. 4 to 6 of this thesis. 

Both Pscaiar and Ptensor are evaluated at around the time of horizon crossing and as such 
they are due to genuine quantum effects. 

The contributions in the third and fourth lines of Eq. (3.5), also dubbed as "classical 
one-loop", can be considered as classical loop contributions arising after the perturba- 
tion modes leave the horizon. The distinction between classical and quantum loops is 



Pc 




(3.7) 



Pscaiar = [9iln{k) + g2] 



(3.8) 
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intended as for example in [109]: quantum loops find their origin in the Lagrangian in- 
teraction terms between the inflaton perturbations and the gravitational modes or from 
self-interaction of (50; classical loops are corrections merely coming from the expansion of 
( using the 6N formula and originate from zeroth order terms in the Schwinger-Keldysh 
formula. 

Finally, the second line of (3.5) includes the integral of B^{ki, /c2, ks), the bispectrum of 
the scalar field defined by 



(50,-;%-^/0,-;) = (27r)35(3)(fc-; + k, + h)B^{k,, k2, ks) 



(3.9) 



and from [13] we have 



V~^HtF{ki) 



(3.10) 



rap 



where mp is the Planck mass, e is the slow- roll parameter (e 
function of the momenta moduli ki of dimension (mass)~^. 



H/H^) and F is a 
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4. Perturbative expansion of the Lagrangian in P(X, 0) theories 



In this and in the next two sections, we will review the computation of the tensor loop 
corrections to Pq. For our purposes, the exponentials in Eq. (2.2) need to be expanded 
up to second order in the interaction Hamiltonian Hi 



T 

if 



0) (4. 



r dri (Ht{ri)-Hj{r^' 

J — oo ^ 

{Q\T[e r dv'(Htiv')-HT{v')) r dri' iHt{rl')-Hj{ri' 

J—OO ^ ^ J —CO ^ 



2 

where 6(t) = 5(j)^^{ri)5(l)^^{ri). One- loop power-spectrum diagrams require an expansion 
of the interaction Hamiltonian to third and fourth order in the field fiuctuations, i.e. 
Hi = nf^ + Hj^\ We provide in Figs. (3) and (4) the diagrammatic representation 
of the leading order corrections that we will find for the diagrams with tensor loops in 
single-field slow-roll inflation. The continuos lines represent scalar propagators, whereas 
the dotted lines indicate tensor propagators. In order to derive this result and the an- 
alytic expressions for these diagrams, we need to first calculate and expand Hi up to 
fourth order in the field perturbations 6(f) and ^7. The starting point is the Lagrangian 
of the system. 

We will begin with a more general Lagrangian for the scalar field than the usual 
= {\/2)g^"dfj^(j)dij(j) — ^(0), by introducing a non-conventional kinetic term, i.e. 

S= - f d^xdty/^ \mIR + 2P(X, 0)1 , (4.2) 



2 

where P = P{X, 0) is a generic function of the scalar field and oi X = ^g^ud^<pd^(p ^ind 
R is the Ricci scalar in four dimensions. Notice that the action (4.2) reduces to the 
standard case ii P = X ~ V, where V is the potential for the scalar field. 
Theories of infiation where the Lagrangian kinetic term is a generic function of the 
scalar field and its first derivatives, like in Eq. (4.2), are string theory-inspired. They 
represent interesting alternatives to the basic infiationary scenario because of their non- 
Gaussianity predictions. The crucial quantity in this sense is represented by the speed 
of sound = (dxP)/ {dxP + 2XdxxP), which is allowed to vary between and 1. 
The perturbative expansion of the interaction Hamiltonian in this kind of models has 
coefficients proportional to inverse powers of the sound speed and therefore, for small 
values of Cg, allows both for non- negligible loop corrections to the power spectrum of 
the curvature fiuctuations [113] and for large values for the amplitudes of three [31] and 
four [41, 42, 43, 44, 45, 46] point functions. In this thesis, we will carry out the calcu- 
lations of the interaction Hamiltonian for these general theories up to a certain point 
and then, for simplicity in the presentation, focus on the canonical case (the remaining 
computations for more general Lagrangians will be found in [113]). 

Let us list the background equations for the system 

2H + 3iJ^ = -P, (4.3) 

12 



3H^ = 2XPx - P, (4.4) 
X {Px + 2XPxx) + 2Vs{2XPx - P f^XP: 



X 

= V2X{P^-2XPx^), (4.5) 
where a dot indicates a derivative w.r.t. cosmic time and, to zeroth order, we have 

The so called flow-parameters are defined as 

e - ~. (4.6) 

These quantities reduce to the slow-roll parameters in the standard case, so it is natural 
to assume |£| <^ 1 and |?7|. It is not correct to talk about slow-roU if P is left as a generic 
function of X and </>, since the smallness of e and rj does not necessarily indicate that 
0^ ^ H'^ and |0| ^ \3H(j)\. It can be convenient to decompose e as the sum e = e^ + ex, 
where 

6 dH 



(4.8) 



-X - (4.9) 

The parameters that are expected to appear in the perturbative expansion of the 
Lagrangian are 

ct = ^ , (4.10) 

' Px + 2XPxx' ^ ' 

CsH 

u=l-\. (4.12) 

J: = XPx + 2X^Pxx, (4.13) 

A = X^Pxx + ^X'Pxxx, (4.14) 

n = X^Pxxx + Ix^Pxxxx, (4.15) 
5 

where Cg is the sound speed. Cg is allowed to vary between and 1, so the quantity |m| 
can freely range between and oo. The only assumption we make is s -C 1, from Cg 
being constant in the standard case. 

4-1. Arnowitt-Deser-Misner (ADM) decomposition for P{X,(f)) theories 

The Lagrangian in Eq. (4.2) will now undergo a perturbative expansion in terms of the 
field fiuctuations 6(f){x, t) = 0(x, t) — 4>o(t) (0o is the homogeneous background value for 
the field) and of the metric fiuctuations. 
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Figure 3. Diagrammatic representation of the leading order (tensor mode) corrections 
from i/j'*' to the power spectrum of 5(p. 



Figure 4. Diagrammatic representation of the leading order (tensor mode) corrections 
from Hj to the power spectrum of Scf). 



It is convenient to adopt the 3+1 Arnowitt-Deser-Misner (ADM) splitting for the metric. 
In the spatially flat gauge the perturbed metric is 

ds^ = -N'^dt^ + h,j{dx' + N'dt){dx^ + N^dt), (4.16) 



h,, = a'{t){e''),j, (4.17) 

where a(t) is the scale factor, is a tensor perturbation with dijij = ■ju = (traceless 
and divergenceless) and det{e'^)ij = 1. Notice that repeated lower indices are summed 
up with a Kronecker delta, so di'jij stands for 6''^di'jkj and 7^ = S'^^'jij. 
In the ADM formalism, the action (4.2) becomes [13] 

S =^J dtd^xVh [nR^'^^ + 2A^P + iV"^ {EijE'^ - E^)] , (4.18) 

where i?*-^-* is the curvature scalar associated with the three dimensional metric hij and 

= \ {hij - \7iN, - \7jNi) , 
E = h'^E.j. 

A dot indicates derivatives w.r.t. time t, all the spatial indices are raised and lowered 
with hij and units of Mpf = 8itG = 1 will be from now on employed. To Ath order we 
have 

R^''^ = -\daaA-fal. (4.19) 

The lapse and shift functions, and N^, can be written as 

A^= 1 + a, 
AT, = dj9 + f3j , 

where a, 9 and /3 are functions of time and space is divergenceless). We have 
exploited the gauge freedom to set two scalar and two vector modes to zero, thus 



14 



leaving one scalar mode from N, one scalar and two vector modes from Nj and two 
tensor modes (the two independent polarizations of the graviton) from hij, together 
with the inflaton field perturbation 6(f). N and Ni are non-dynamical degrees of freedom 
and can be expressed in terms of the other modes (50 and 7jj), once the Hamiltonian 
and the momentum constraints (we derive them in the next section) are solved. 



4.I.I. Solving Hamiltonian and momentum constraint equations 

Momentum and Hamiltonian constraints are derived from varying the action w.r.t. the 
shift and lapse functions respectively. It turns out that, in order to expand the action 
to a given order n, it is only necessary to perturb N and Ni up to order n — 2 [13, 31]. 
Therefore we will solve the constraints to second order in the metric and scalar field 
fluctuations. 

Let us begin with the expansions 



a 

e -- 



■■ «1 + 02, 
01 + 02. 



where ai and 02 are respectively first and second order in the fields fluctuations 
(similarly for Pn and /32i, and for 61 and 62). Let us then expand P to second order. P 
is a generic function of X and 0. We first need the expansion of X 



X 



g'^d,^ 

1 r 
2 



Ar-2 

= Xo + AX 
where X* 



2 

+ 50 



+ 2c/°^ai00 + g'^d, 
' + 2N-'^N'di5(t) (0 + 50 



X2 



di6(j)dj 



(4.20) 



: W^Nj, Xq is the zeroth order part, i.e. Xq = ^ and AX is the 
perturbation to the desired order (AX = AXi + AX2 + AX3 + ...). Notice that 
0(t, x) = 0o(i) + 50(t, x), but for simplicity we will suppress the subscript '0' in the 
background value of the field. 
The expressions for the perturbations AXj become 



AXi = 2X0 



AX2 = Xo 



—. «! 



- 4ai4- - 2a2 + ^ai^ - 2N[d~ 



-' 0^02 



di5(j)di6(j) 



and so on for AX3 and higher order terms. The expansion of P{X, 
order becomes 



(4.21) 

(4.22) 
up to second 



P(X, 



1 



1 



Po + PxloAX + P^|oA0 + -Pxxlo(AX)^ + 7^P00|oA0' + Px^|oAXA0 



21 



(4.23) 

where as usual the subscript '0' indicated the zeroth order, Px = dxP, P4) = d^P and 
similarly for the second order derivatives, A0 = 50 and AX needs to be expanded up 



15 



to the needed order. 



We are now ready to write the momentum and Hamiltonian contraints 



(4.24) 
(4.25) 

The momentum constraint to first order reads 

2Hdjai - ^d^Pi, = Px4>dM , (4.26) 

where H = a/a is the Hubble parameter. Eq. (4.26) can be solved to derive ai. Taking 
the derivative of both sides of (4.26) and using the divergenceless condition for /3, we 
have 



2H 



(4.27) 



Using the solution found for ai, we find d'^Pij = 0, from which we can set = 0. Here 
^2 ^ S^^didj, which we will indicate in the rest of the thesis also as didi, and from now 
on we define /3j = /32i for simplicity. 



The momentum constraint to second order is 

2Hdja2 - AHaidjai - ^djUid'^Oi + -^diaididjOi - \diai'ii^ 

a^ a^ ■'2 

+ -^iikdakj - -^likdiikj - -^iikdjjik + -^diOid^-fij 



Pxd, 



2XPxxd, 



2XPxx4"^idj6(f) — Px4>0Lidj 



+ 



50. 



The solutions are 



a; 



012 = — 

2 



2Ha'^ 



1 



Aa'^H 



+ [djiikdj-fik 

XXg.2 



XP 



where S = 9' 



H 

dj6(f)dj64>, and 



2H 

+ dj6(pdjC 



didjOid'^'jij 
- {djaidjScf) + 



1 



-^d'^(5j = 2Hdja2 — AHaidjUi -djaid'^Oi + -^diOididjOi 



2a 



a" 

ld.a,i., + \i.Aik, - \i.kdak, - \i.kd,i.k 
1 

32 



^diOid jij - Pxdj 
2XPxx^aid,5<P + Px 
16 



2XPxxdjS(p5(j) 



(4.28) 



(4.29) 



(4.30) 



Let us now move to the Hamiltonian constraint which provides 



AH 



'5<p 



X 



'5<p 



(4.31) 



to first order and 



4H 



^ 8^02 = (-2ai) kxPx^ + 20Pxx^'^ + 2XPx<^50 + 8PxxxX^^-i 



AH 

+APxx^XH(t>+—d^O, 

(X 



a?(t) a* V 



%d,dA + -{d,dje^Y +[-6H' + 2XPx + 4X'Pxx) 3«f - 2^2 

6<p 



+4al {sX^Pxx + 2X^Pxxx) + [2XPx + IGX^Pxx + SX^Pxxx 



AXPxcj) + SX'^Pxxij) 



AX'Pxx - 2XP 



X 



+2XPx(j,(j, 
to second order. 



AH 



iijiij + -^da'jigda'jig YlijdidjOl 



(4.32) 



4.1.2. Reduction to the canonical case 

In the canonical case, to zeroth order in perturbation theory P = (0^/2) — V{(j)), so 
Px = 1 and P^n = —d"'V/d(f)"' with all other derivatives of P being zero. The solutions 
above therefore reduce to [56] 
1 • 



2H 



AH 

a 



2 d% = -2V^6<j) 



2ai (-6H^ + (0)^ 



(4.33) 
(4.34) 



"1 1 o-2v^ 1 



1 



8H 



(4.35) 



AH 



92^2 = 2ai 



■ ■ AH 

2050 + —d^e 



a'- 



'^j^mdiScj) - -didjOididjO^ (4.36) 



a' 



3«? - 2«2) (0^ - 6/7^) - (50^ - — 9i(509i50 



4a2 



^d^/3j = \d'^aidjd'^ei - \d,ndjaidmd^ei + \dmaid^djd'^ei 
2a^ 
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- ^djaid'^Oi - \dmdjdiaididmOi + ^d'^diaididjOi 

1 1 

-didjaidid^6i + —drrAaidmdidj9i + dmdj6(pdrn6(f) 

— d'^6(j)dj6(f) + djS4>d'^6(p — dm8(t)dmdj54> — -d^ (iijdiai) 

+ {iikdau,) . (4.37) 

where V^^ = d'^V/dcjP' and is the inverse of the laplacian operator. Notice that the 
equations (4.33) through (4.37) agree with the resuhs obtained in [152], if we set to 
zero. 



4.2. Fourth-order expansion of the Lagrangian in the canonical case 



The expansion of the action up to 4th order can be now derived by plugging the solutions 
(4.27) through (4.32) in (4.18). The final expression is quite lengthy and can be found 
in [113]. We will here only report the 4th order expansion of the action in the canonical 
case 



^4 



dtdr'xa^ 



1 

24 



—dj9idj5(l)dm.9idr, 



+ 



a' 



a\a2 



-OLr 



2V^ai64> + ai - 



2 



dj92 + Pj)dj5(j) 



2a2 



^Iiklkjdj6(j)di6(l) - ai'yijdj6(f)di6(l) + a2di6(f)di6(f) 

- (pdj6(j) {jijdi92 + jij(3i + jijdi9i) - dk92iabdblak - hiabdblak 

- -iabPkdklab - ai[H-1abdadb92 + iabdadb92 + labdaPb) + ]pabdkl abdk92 

+ ^( - 8'yipdidj9idpdj92 - 4'jipdidj9idpl3j - 4'yipdpdj9idjPi 

- dq9idajgdidj92 - dq92dajgdidj9i - /3gdajqdidj9i 

- dq9idajgdil3j - dg9idajgdjl3i + dg9idg'yijdidj92 

+ 2dg92dg-fijdidj9i + 2dg9idg'yijdi/3j + 2/3gdg-fijdidj9i 



(4.38) 



Similarly, the loop computation will be from now on performed considering this simpler 
case. 
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5. Tensor loop corrections to 



Let us then consider tlie terms in the interaction Hamiltonian Hj that involve the tensor 
modes. The third order action in single-field slow-roll infiation with non-zero graviton 
fiuctuations was calculated in [13]; we will focus on the leading order term in the slow-roll 
parameters, so we have 



d x^ijdi5(j)dj5<j). 



(5.1) 



The fourth order action is given by Eq. (15.2.1). Notice that some of the interaction 
terms involving the tensor modes in (15.2.1) appear with time derivatives, therefore the 
construction of the path integral formula requires additional care compared to the case 
where time derivatives only appear in the kinetic term of the Lagrangian. This issue will 
be discussed in Appendix 15.1. Also, it is possible to show that in Eq. (15.2.1), of all the 
leading terms in the slow-roll expansion, only one will provide a non-zero contribution 
to the loop correction (see Appendix 15.2 for a detailed analysis) and contribute to the 
interaction Hamiltonian to fourth order which becomes 



4 



(5.2) 



where the tensor fiuctuations are 
with 



A,A' 



(5.4) 



The equation of motion for the eigenf unctions 6(j)k{t) can be derived in the approximation 
of de-Sitter space from the second-order action 



drj 



(5.5) 



(where drj = dt/a{t) is the conformal time) and they are given by the well-known 
expression 

H , 



Ukiv) 



'2k^ 



ikr]) e 



-ikr) 



(5.6) 



In the same approximation, the eigenfunctions for the tensor modes '^kiv) given by 
ul = 2uk. 



Let us now begin with the one-loop one-vertex part of the diagram (given in Fig. 3) 
which we label with the subscript (IL, Iv); this can be written as [104], [105] 



{s<PkSv*)s<f>kSv*))( 



lL,lv) 



drj 



H 



(4) 



(V),50,q(r/*)50,-;(r/*) 



(5.7) 
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We will study this in detail 



X 5(3) (^-; + g-)5{3)(/r + + q^) + c.C, (5.8) 



where the extra factor of 2 accounts for the number of equivalent diagrams obtained by 
permuting the field contractions, Ui^ij]) is given by Eq. (5.6) and 

P^j = E (^) = 2 sin^ 95,,. (5.9) 

a,a' 

Integration and the use of the delta function lead to a simpler form 

X (1 — ikrj ) (1 + iqrj )(1 — iqrj ) + c.c, (5.10) 

This equation is exact except for the approximation of using the de Sitter space for- 
mula for the scale factor, 0(77) = — (ifr/)"^, and evaluating the Hubble radius H{r]) at 
the time 77*. The reason why this is allowed is the following: the contribution to the 
integral w.r.t. time from regions well before horizon crossing is negligible compared to 
the contribution due to the region around horizon crossing [13, 104, 105]; in addition to 
that, we are choosing 77 to be just a few e-folds after horizon crossing, so we can assume 
that the Hubble radius (as well as any of the slow-roll parameters of the theory) will 
not undergo a big variation during this interval of time. The same approximation will 
be applied to the two-vertex diagrams. 

We first solve the time integral. It is convenient to perform a change of variale like 
in [106], i.e. we set x = —krj' and x* = —krj* so that 

= S^'\k, + r,)g / ^sin^^ Im[£ ^^e-(^'--) 

X (1 + tx*?(l - tx'fn + - 4^')] ■ (5.11) 

k k J 

After integrating the imaginary part, we end up with the following result 

Ht fd^Q. 2^2A;2(3 + x*2) + g2(5 + 5a,*2^2x*4) 



= ^^'H^"; + ^"^2)^27r^[^(3 + X*') / y + ^(5 + 5x*' + 2x*') J dqq 

where the factor 4/3 comes from integrating with respect to the azimuthal angle 6 (notice 
that that the reference frame in momentum space has been chosen in such a way that 
the external wave vector k lies along the positive z axis). We now solve the momentum 
integrals. Both the logarithmic and the quadratic one exhibit ultraviolet divergences 
and the logarithmic part diverges also at very low momenta. Ultraviolet divergences 
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(5.12) 



can be treated as in flat space; the infrared logarithmic divergence is flxed introducing 
a momentum lower cutoff to be interpreted as a 'box size' [110, 111, 117, 118, 119] 
which can be flxed to be not much larger than the present horizon [110, 111]. As an 
example, consider the flrst integral of Eq. (5.12) which is convenient to split as follows 
dq dq dq 

if-i q i^-i q Jk q ^ 
where we have introduced an upper cutoff A. The flrst integral gives ln{k£); the 
second integral can be renormalized introducing a counterterm —In^A/ko), where feg 
is a renormalization constant. The flnal result for Eq.(5.12) reads 

{HkM)Hk,{v*))iiLM = 7r6^'Hh + k2)'^{S + x*')[\n{ki) - \n{k) + a 
where a is a left over constant from renormalization. 



Let us now focus on the one-loop contribution from the 3rd order action with the 
gravitons (see Fig. 4 for its diagrammatic representation) 



X 



T 



dv" (Hnv")-HY{v" 



T 



6<P,^^{v*)S<P^^{r]*){A + B + C + D 



where 



n* , . rv* „ . 
A= driHl dr] Hf , 

J —oo J—oo 

rv* , rv' I, 
B= dr]Hj dr] Hf, 

J — oo J—oo 

rv* , . rv* n 
dr] Hj / di] HJ 



-oo 
V* 



oo 
V* 



c = 

D = - I dr]Hj I dr] H+ 

J—oo J—oo 

It is easy to check that B = A* and C = C* = D. We can write Eq. (5.13) as 

dr] fv* dr]" 



X 



(5.13) 

(5.14) 
(5.15) 
(5.16) 
(5.17) 

(5.18) 



{Hr]'y J-oo {Hr]" 

where the factor sin^ 9 comes from contractions of the polarization tensors with external 
momenta [120] 

/ ^ ^\ 2' 

q- k 



eij{q)k'k^ 



7i 



1 - 



qk 



V2 



sin" 9, 



(5.19) 
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and the wave fuctions Wf are 

wf{v',v") = Uk{vlul{ri)Uf^{r]*)ul{rl')[u^^_^{ri)u*^^^ 

xQiji - r]") + u*f:_^{r]')u^^_^{r]")u*{T]')Ug{r]")e{r]" - r/')] , 

X Q{r]" -ri)+ u*^_^{ri)u^^_^{ri')ul{ri)Uq{ri')Q{ri - V')] , 
w'j^iv\v") = -%^(^*)Mfc(V)Mfc(r/*)Mfc(V')Mffc_g-i(V)M|g_g-^(V')<(V)Mg(V') , 

so w^{-r]\vi") = wf{rj'\rj') and is a real number and 'wf{ri\rj") = wf*{7]" ,7]'). We will 
label the two contributions by A and C, so that the one-loop contribution with two 
vertices to the two point function will be broken into two parts 

(50,;;(r/*)50,;;(r/*))(i^,2.) = (%-l(^*)50.:;(^*))fiL,2.) + (%:l(^*)^0fc;(^*))fiL»- (5-20) 
Let's look in details at the two parts. 

{HuS^*)5^kM))tiLM = -S^'\ki + k2)^2 I ^^^,e-'^''\l + ^kv*f (5.21) 

ZK J q \k-q\ 

X r %e''''^''-'^-^''-'^\l-ikr]'){l + iqri'){l + i\k-q\ri') 

J-oo 7]' 

X r %e*('=+''+l^^-'^^)''"(l-2A;r/")(l-igV')(l + ¥-^V') + c.c. 
J-oo rj" 

The second time integral has e*^'' — 4?- + ^77" — i (^p^) its primitive function, where 

g = k + q+\k — q\,h = —qk — {q + k)\k — q\ and c = qk\k — q\. This should be evaluated 
between —00 and 77'. It is soon evident that the lower bound represents a problem for 
this evaluation. We need to remind ourself, though, that the choice of the integration 
time contour needs to be deformed and to cross the complex plane to account for the 
right choice of the vacuum [13]. This is done by integrating in a slightly imaginary 
direction, i.e. taking rj" — )■ r]" + ie|?7"|, where e is a fixed small real number; so for 
example 

/ dr/ e'^''' = V- (5-22) 
J-oo ik 

With this contour prescription, our integral in 77" vanishes at —00. Performing the same 
change of variables as in (8.14) 

y |/c — g| 

( k c t f gb — c\\ f d I s ,2 

(5.23) 
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IX 



where c? = q—k+\k—q\ and s = kq+{k—q)\k—q\. The result of the integration w.r.t. time 
is a polynomial function of sin2x*, cos 2a;*, Si(2x*), Ci(2x*) and their products with 
coefficients which depend on b, c, d, s, k. Notice that in the large scale limit x* — )■ 
a singularity similar to the one found in [106] shows up in our result (see also [121] for 
a recent discussion on these kind of singularities). However, by evaluating the power 
spectrum of 6(j) just a few e-folds after horizon crossing, we are safe from these kind of 
singular behaviour [109]. 

The next step consist in performing the momentum integral. The integrals we need to 
evaluate are of the following kind 

d^q sin^O ^ , . „ 
-T— -73/(?1> (5.24) 
y \q — k\ 

where /(g) is a sum of functions of momentum. Let us begin for simplicity by considering 
the constant term of the sum, i.e. let us study 

d^q 



\q-kf' 



(5.25) 



For the specific case of equation (5.25) the integrand function has singularities at g = 
and at q = k and shows no ultraviolet singularities. Based on an approximate evaluation 
performed considering a sphere of radius around q = 0, where <^ k, the integral 
is proportional to a function ln{ki). The same result can be obtained working in a small 
sphere around q = k after a change of variables go = (1~ k. The contribution from large 
values of q is negligible w.r.t. the ones from the singular points, so the integral over the 
whole momentum space is expected to be proportional to ln{ki). The exact value of 
the integral can been found after a change of variable from the (g, 6) to the {q,p) space, 
where p = |g - k\ and is equal to (167r/225P) (1 + 301n(H)) ~ k"^ (10^^ + lOln(H)). 

Integrating Eq. (5.23), we find ultaviolet power law and logaritmic singularities in 
addition to infrared logaritmic contributions. The final result of the integration is a 
function of x* = e~^*, where A^^, = ln(a*/afc) is the number of e-foldings from horizon 
crossing 

(50fc;(^*)%-;(^*))fiL,2.) = + k2)^{a, \n{k) + a,\n{ki) + a-,) , (5.26) 



where ai, 02 and as are functions of x* (see Appendix 15.3). We are calculating the 
two point function for the scalar field a few e-foldings after horizon crossing, so x* may 
be chosen to range between 10~^ and 10~^. In this range Oi ~ C>{1) and negative, 
as = -16/(15x*2) + (8/15)(5 - 8Ci(2x*)) and ag = -8/(225x2) + C(l) + p, where p is 
a left-over scheme- dependent renormalization constant of the kind present in equation 
(5.13). 
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Let us now move to part C of Eq. (5.13) which we give below 



X 



-oo 7] 



V* drj" 



-igr] 



oo r] 

3^ „;„4z 



1 + {kr]^ 



Q*iv" 



*\2 



Ht fd-'q sin^e 



\k-qf 



1 + {kri' 



*\2 



X 



i?e / drj e 



tgr) 



Im / dr] e 



^gr, Qijl 



(5.27) 



where Qiji) = 1 + igrj' + brj''^ — icr]'^ . Let us integrate over conformal time 



°° dx k"^ 

72 



X* 



k X 



• 9 



g I h ,2 c /3 
k k^ k-^ 



e 



fc2 



fc^ cfc ^ (gb — c)k'^' 
\ X +1 



X 



9 



9^ 



(5.28) 



where again the integration has been performed by continuing rj' to the complex plane, 
i.e. {rj' T]' + ie\r]'\), and then taking the limit e — )■ 0. 
We are now ready to integrate over momentum 



Ht fd^q sin'^e 



2P 



l' \k-q\' 



1 + x* 



k' 



2k^c 2k%c k^c^ k%^ 



„*2 



+ 



+ 



(5.29) 

2;.. g g.. g^ g^ g^ ^ 

Similarly to what we have done in part A, one can check that there are no ultraviolet 
singularities in the remaining five integrals although some infrared logarithmic 
contributions are still present and the final result is 



(5.30) 



where ci = (1/225) (8/a;*^ + 107 + 50^*^) and ca = (16/15x*2) + (4/15). Notice that 
the (x*)~^ coefficients in ci and C2 exactly cancels the (a;*)~^ coefficients in 02 and 03. 
This is not surprising: based on [106, 109], we expect we might observe a logarithmic 
singularity if we push x* — )■ in our results (which is indeed present in the Ci(2x*) term 
of 02), but no power-law singularities are actually expected. 
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6. Complete expression for at one loop 

Let us now collect our results in the final formula for the power spectrum of the curvature 
perturbation ( computed up to one-loop level. This can be derived from Eq. (3.4), which 
follows from the 6N formula. Summing the main results of the previous section, Eqs. 
(5.13), (5.26) and (5.30), we obtain the one-loop graviton correction to the infiaton 
power spectrum 

{6<P^S'^*)S<l>,^^iv*))iL = 7T6'~'\h + h)^ [/i ln(A;) + f, \n{ki) + fs] , (6.1) 
where 

/i = -^(25 + 15x*2 + 4x*4) , (6.2) 

/2 = 2 + {2/3)x*^ + 02 + C2 , (6.3) 

and /s is given by a left-over scheme-dependent renormalization constant plus 
contributions of order 0{1) (see Appendix 15.3 for the complete expressions of a2, 
C2 and /s). If we calculate the two point funtion of 6(f) a few e- foldings after horizon 
crossing, i.e. x* ranges for example between 10"^ and 10"^, /i reduces to a negative 
constant of order 0(1) and /2 ~ 4 (1 — Ci(2a;*)) ~ C(10). In the limit where a;* — )■ 1 
both fi and /2 turn out to be of order unity. 

In order to understand which is the dominant contribution in Eq. (3.5) and how big 
it is, one needs to (i) know the slow-roll order of the coefficients A^*^*); A^*^^) ~ e^^^^, 
NC^) ~ A^(3) gi/2. compute the integrals involving the power spetrum P{q). 
This is discussed in details in Ref. [109] (for the case of scalar perturbations only), see 
in particular Sec IV of [109]. It turns out that the crucial quantity is represented by 
the number of e-foldings of inflation between the times of horizon exit of the mode 
which corresponds to the infrared cutoff, and the time of horizon exit of the mode k we 
want to observe. However, to deal with observable quantitites one has to choose i not 
much bigger than the present cosmological horizon Hq^ [HO, 111]. 
The relevant point about Eq. (6.1) is that it gives in Eq. (3.5) a contribution which is 
of the same order of magnitude as those coming from loops which accounts for scalar 
perturbations only. Since in terms of the slow- roll parameters (^N^^^^ ~ e*~^ the 
magnitude of the one-loop graviton correction turns out to be 

APl'°°^{k)r^^a{k)^V:{k), (6.4) 

where we have used Eq. (3.6) for the power spectrum of the inflaton fleld. In Eq. (6.4) 
a{k) includes the various coefficients of Eq. (6.1), and it is 0{1). Eq. (6.4) allows a 
more direct comparison with the results of Ref. [109], showing that the graviton con- 
tributions to the one-loop corrections are comparable to the ones computed only from 
scalar interactions. Notice that also for the tensor contributions we find terms of the 
form In(fc). 
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Summarizing our work: beyond linear order, the tensor perturbation modes produced 
during inflation unavoidably mix with scalar modes; this fact alone would require to 
include the tensor modes for a self-consistent computation. Most importantly, despite a 
naive expectation suggested by the fact that the power spectrum of the tensor modes is 
suppressed on large scales with respect to that of the curvature (scalar) perturbations, 
our results show explicitly that their inclusion is necessary since their contribution is not 
at all negligible with respect to the loop corrections arising from interactions involving 
the infiaton field only. 
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7. Study of perturbations in anisotropic cosmologies 

The Bianchi models represent a classification of all homogeneous and anisotropic 
cosmologies. In the cosmic no-hair conjecture, any initial background, in the presence 
of a positive cosmological constant, eventually evolves in a de Sitter universe, where all 
the initial existing anisotropics are rapidly washed out. This conjecture was proven to 
be true for all of the Bianchi models except for the Bianchi-lX by Wald in [72]. For the 
Bianchi-IX this result is also true under the assumption that the cosmological constant 
overcomes the spatial curvature terms. More recently, there has been a revived interest 
in homogeneous but anisotropic models of the early Universe and several attempts to 
understand what kind of observational signatures these models might produce. A close 
look has in fact been given to the fluctuations in the metric and in the energy tensor 
developing during a hypotetical anisotropic stage either in a pre-infiationary epoch 
[78, 79] or during inflation itself [76, 77, 144]. The Bianchi-1 model is the simplest 
of all anisotropic and homogenous models and has often been the choice for such kind 
of studies. In the Bianchi-I model, the metric has a form 

rfs^ = _ / 2/3\ ^^i^^j ^^^^^ 

V / ij 

where a is a function of time as well as Pij. The latter is a 3 x 3 diagonal traceless 
matrix that anisotropizes the volume expansion 

{e'^)^. = 6,,e'^% EA = 0. (7.2) 

i=l 

A metric described by Eq. (7.1) but with constant functions a and (3 is often presented 
in the form 

ds^ = dt^ - t^P'dxl + t^v^dxl + t^'P'dxl, (7.3) 
where the coefficients pi, p2, Pa are numerical constants satisfying the relation 

Pi +P2 +P3 =Pl +P2 +P3 = 1- (7-4) 
The three parameters can be equal in pairs in the cases ( — |, |, |, ) and (0, 0, 1). In all 
other cases they are distinct, one being negative and the other two being positive. The 
universe described by (7.3) is spatially flat and, for any possible value of the coefficients 
Pi, its volume element is equal to yj—g'-^'d^x = td^x, where is the determinant of the 
three- metric. 

In [78], we consider a pre-inflationary epoch characterized by an initially expanding 
type-I Bianchi universe that evolves with an energy density not yet dominated by a 
cosmological constant. If by the time the cosmological constant has become the domi- 
nant form of energy the metric is still homogeneous, the Universe will eventually enter 
a de Sitter epoch; the fluctuations of this metric and their growth, however, can play 
an important role in determining the evolution of the Universe at this early stage. The 
question we ask concerns the kind of initial conditions that are needed for the Universe 
to remain homogeneous up to the time the cosmological constant eventually dominates 
the total energy density. We employ a metric as in Eq. (7.1) for the background and 
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solve the first order Einstein equations for tfie perturbations considering a pressureless 
fluid. As expected, the anisotropy in the background is responsible for a coupling at the 
linear order of one of the tensor modes with the density contrast, so there is a correlation 
between the scalar and the gravitational perturbations. Moreover, the evolutions of the 
two tensor modes, the "free" one and the one that is coupled to the scalar, are very 
different from each other and very much depend on the scale of interest. It turns out 
that, for a reasonably large set of initial conditions, the growth of these perturbations 
could be fast enough so as to lead the Universe into an inhomogeneous state before the 
cosmological constant becomes the dominant form of energy. 

A similar scenario was investigated in [79], where both expanding and contracting 
Bianchi-I Universes are studied during a pre-infiationary era. A scalar field dominates 
the matter content. In particular, the evolution of the gravitational perturbations of 
the metric (not to be confused with the gravitational waves generated from the vacuum 
fluctuation during inflation) is analysed, to point out that they may grow to significantly 
alter the geometry of spacetime before inflation begins and, in any case, they may leave 
an imprint on the observed CMB power spectrum. 

Also in [76, 77] a perturbative analysis was performed for Bianchi-I cosmologies lead- 
ing to anisotropic inflation; in addition to that, predictions for the power spectrum of 
curvature fluctuations and gravity waves produced during inflation are provided for this 
models. Finally, a "vector hair" model and its observational consequences in terms of 
anisotropic signatures in the cosmological correlation functions and primordial gravita- 
tional waves, were recently discussed in [144]. 
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8. Inflation and primordial vector flelds 

The attempt to explain some of the CMB "anomalous" features as the indication of a 
break of statistical isotropy is the main reason behind ours and many of the existing 
inflationary models populated by vector fields, but not the only one. The first one of 
these models [122] was formulated with the goal of producing inflation by the action of 
vector fields, without having to invoke the existence of a scalar field. The same motiva- 
tions inspired the works that followed [123, 124, 125]. Lately, models where primordial 
vector fields can leave an imprint on the CMB have been formulated as an alternative to 
the basic inflationary scenario, in the search for interesting non-Gaussianity predictions 
[90, 91, 92, 93, 94, 95, 96, 97, 98]. Finally, vector fields models of dark energy have 
been proposed [126, 127, 128, 129, 130, 131]. All this appears to us as a rich bag of 
motivations for investigating these scenarios. 

Before we quickly sketch some of them and list the results so far achieved in this di- 
rection, it is important to briefly indicate and explain the main issues and difficulties 
that these models have been facing. We will also shortly discuss the mechanisms of 
production of the curvature fluctuations in these models. 

Building a model where primordial vector fields can drive inflation and/or produce the 
observed spectrum of large scale fluctuations requires a more complex Lagrangian than 
the basic gauge invariant C^ector = "{yj—g /^)F^" F^iu- In fact, for a conformally in- 
variant theory as the one described by H vector-, vector fields fluctuations are not excited 
on superhorizon scales. It is then necessary to modify the Lagrangian. For some of 
the existing models, these modifications have been done to the expense of destabilizing 
the theory, by "switching on" unphysical degrees of freedom. This was pointed out in 
[132, 133, 134], where a large variety of vector field models was analyzed in which lon- 
gitudinal polarization modes exist that are endowed with negative squared masses (the 
"wrong" signs of the masses are imposed for the theory to satisfy the constraints that 
allow a suitable background evolution). It turnes out that, in a range of interest of the 
theory, these fields acquire negative total energy, i.e. behave like "ghosts" , the presence 
of which is known to be responsible for an unstable vacuum. A related problem for some 
of these theories is represented by the existence of instabilities affecting the equations 
of motion of the ghost fields [132, 133, 134]. 

In the remaing part of this section, we are going to present some of these models to- 
gether with some recent attempts to overcome their limits. 

In all of the models we will consider, primordial vector fields fiuctuations end up either 
being entirely responsible for or only partially contributing to the curvature fiuctuations 
at late times. This can happen through different mechanisms. If the vector fields affects 
the universe expansion during infiation, its contribution to the total can be derived 
from combining the definition of the number of e-foldings (A^ = / Hd£) with the Einstein 
equation [H^ = (87rG/3)(p0 + pa), Pa being the energy density of the vector field and 
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the inflaton energy density) and using the (5N expansion of the curvature fluctuation 
in terms of both the inflaton and the vector fields fluctuations (see Sec. 9). To lowest 
order we have [92] 

U = ^5 A, (8.1) 

Imp 

where a single vector field has been taken into account for simplicity (mp is the reduced 
Planck mass, A is the background value of the field and 5 A its perturbation). When 
calculating the amplitude of non-Gaussianity in Sec. 10, we will refer to this case as 
"vector inflation" for simplicity. 

A different fluctuation production process is the curvaton mechanism which was initially 
formulated for scalar theories but it is also applicable to vectors [135, 136]. Specifically, 
inflation is driven by a scalar field, whereas the curvaton field(s) (now played by the 
vectors), has a very small (compared to the Hubble rate) mass during inflation. Towards 
the end of the inflationary epoch, the Hubble rate value starts decreasing until it equates 
the vector mass; when this eventually happens, the curvaton begins to oscillate and it 
will then dissipate its energy into radiation. The curvaton becomes responsible for a 
fraction of the total curvature fluctuation that is proportional to a parameter, r, related 
to the ratio between the curvaton energy density and the total energy density of the 
universe at the epoch of the curvaton decay [92] 

a = (8.2) 

3 Pa 

where r = 3pA/(3pyi + 4p(^). Anisotropy bounds on the power spectrum favour small 
values of r. 

From Eqs. (8.1) and (8.2) we can see that, dependending on which one of these two 
mechanisms of production of the curvature fluctuations is considered, different coeffi- 
cients will result in the expansion (see Eq. (9.4)). 

In this section we will describe both models where inflation is intended to be vector-fleld 
driven and those models in which, instead, the role of the inflaton is played by a scalar 
field, whereas the energy of the vector is a subdominant contribution to the total energy 
density of the universe during the entire infiationary phase. 



8.1. Self-coupled vector field models 

A pioneer work on vector field driven infiation was formulated by L. H. Ford [122], who 
considered a single self-coupled field A^ with a Lagrangian 

K,,t,r = --^F^,F^"' + V{^) (8.3) 

where F^^y = d^j^By — dyB^j, and the potential ^ is a function oiip = BaB". Different sce- 
narios of expansion are analyzed by the author for different functions V. The universe 
expands anisotropically at the end of the inflationary era and this anisotropy either 
survives until late times or is damped out depending on the shape and the location of 
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the minima of the potential. 



The study of perturbations in a similar model was proposed by Dimopoulos in [135] 
where he showed that for a Lagrangian 

Sector = -If^uF^-" + ^m'B^Bf^ (8.4) 

and for ^ —2H^, the transverse mode of the vector field is governed by the same 
equation of motion as a light scalar field in a de Sitter stage. A suitable superhorizon 
power spectrum of fluctuations could therefore arise. In order to prevent production of 
large scale anisotropy, in this model the vector field plays the role of the curvaton while 
inflation is driven by a scalar field. 

8.2. Vector-field coupled to gravity 

The Lagrangian in Eq. (8.4) may be also intended, at least during inflation, as including 
a non-minimal coupling of the vector field to gravity; indeed the mass term can be 
rewritten as 

L,ector^l{ml + ^R)B^B'^ (8.5) 

where, for the whole duration of the inflationary era, the bare mass mo is assumed to 

r / . \ 2" 

be much smaller than the Hubble rate and the Ricci scalar R = —6 



approximated as i? ^ — 12iJ^. For the specific value C, = 1/6, Eq. (8. 



7 + 



can be 
is retrieved. 



For the Lagrangian just presented, Golovnev et al [123] proved that the problem of 
excessive anisotropy production in the case where inflation is driven by vector fields can 
be avoided if either a triplet of mutually orthogonal or a large number N of randomly 
oriented vector fields is considered. 



The Lagrangian (8.5) with ^ = 1/6 was also employed in [136], where inflation is scalar- 
field-driven and a primordial vector field affects large-scale curvature fluctuations and, 
similarly, in [137], which includes a study of the backreaction of the vector field on the 
dynamics of expansion, by introducing a Bianchi type-I metric. 

8.3. Ackerman-Carroll-Wise (ACW) model 

A model was proposed in [138] where Lagrange multipliers (A) are employed to determine 
a fixed norm primordial vector field B^B^ = 

Lvector D A (5^5^ - m^) - PA (8.6) 

where pa is a vacuum energy. The expansion rate in this scenario is anisotropic: if we 
orient the a;-axis of the spatial frame along the direction determined by the vector field, 
we find two different Hubble rates: along the x-direction it is equal to 

Hi = -^^r^y (8-7) 
nip P{fi) 

31 



and it is given by Ha = (1 + cfi'^)Hf, along the orthogonal directions; fi = m/mp, P is 
a polynomial function of fi and c is a parameter appearing in the kinetic part of the 
Lagrangian that we omitted in (8.6) (see [138] for its complete expression). As expected, 
an isotropic expansion is recovered if the vev of the vector field is set to zero. 

8.4- Models with varying gauge coupling 

Most of the models mentioned so far successfully solve the problem of attaining a slow- 
roll regime for the vector-fields without imposing too many restrictions on the param- 
eters of the theory and of avoiding excessive production of anisotropy at late times. 
None of them though escapes those instabilities related to the negative energy of the 
longitudinal modes (although a study of the instabilities for fixed-norm field models was 
done in [139] where some stable cases with non-canonical kinetic terms were found). As 
discussed in [132, 133, 134], in the self-coupled model a ghost appears at small (com- 
pared to the horizon) wavelengths; in the non-minimally coupled and in the fixed-norm 
cases instead the instability concerns the region around horizon crossing. 

Models with varying gauge coupling can overcome the problem of instabilities and have 
recently attracted quite some attention. In [90], the authors consider a model of hybrid 
inflation [140, 141, 142, 143] with the introduction of a massless vector field 

L D 1 {d,<pd'^<P + d.xd'^x) - \f\4>)F,uF,. + V{<P, X, B^) (8.8) 

where (f) is the infiaton and x is the so-called "waterfall" field. The potential V is chosen 
in such a way as to preserve gauge invariance; this way the longitudinal mode disappears 
and instabilities are avoided. 

Similarly, Kanno et al [144] consider a vector field Lagrangian of the type 

Lvector = -^f{<P)F>^''F,,, (8.9) 

but in a basic scalar field driven inflation model. Very recently, in [145, 146] the linear 

perturbations in these kind of models have been investigated. 

Finally, in [92, 96] varying mass vector field models have been introduced 

Lvector = -\f{ct>)F^^''F^, + ^m^E^E^ (8.10) 

where / ~ a" and m ~ a (a is the scale factor and a is a numerical coefficient). The 
special cases a = 1 and a = —2 are of special interest. In fact, introducing the fields 
Af^ and A^, related to one another by A^ = fB^ = aA^ {A^ and A^ are respectively 
the comoving and the physical vectors), it is possible to verify that the physical gauge 
fields are governed by the same equations of motion as a light scalar field in a de Sitter 
background. Vector fields in this theory can then generate the observed (almost) scale 
invariant primordial power spectrum. 
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8.5. SU(2) vector model 



Let us consider some models where inflation is driven by a scalar field in the presence 
of an SU(2) vector multiplet [97, 98]. A fairly general Lagrangian can be the following 



S = J d'^Xy/^ 



TTlpR 



gtJ-agUl3 



E 

a=l,2,3 



rpa rpa 



M2 ^ 

^ a=l,2,3 



fS.lll 



where is the Lagrangian of the scalar field and F^j, = dp^B^ — d^B^^ + QcE'^^^'B^^B'^ 
{qc is the SU{2) gauge coupling). Both / and the effective mass M can be viewed as 
generic functions of time. The fields 5^ are comoving and related to the physical fields 
by = (i?Q, B^/a). The free field operators can be Fourier expanded in their creation 
and annihilation operators 



5A^{x,v) 



d^q 



E 



e,^(g)af M",(g,r/) + ef (-g) (a%^) ' Mr (g, r/)J ,(8.12) 

\=L,R,long 

where the polarization index A runs over left (L), right (i?) and longitudinal {long) 
modes and 



a, A 



(8.13) 



Here r) the conformal time {drj = dt/a(t)). Once the functional forms of / and M have 
been specified, the equations of motion for the vector bosons can be written. We compute 
cosmological correlation functions up to fourth order considering an action as in (8.11). 
The expression of the correlators that we derive, prior to explicitating the wavefunction 
for the gauge bosons, apply to any SU{2) theory with an action as in (8.11), both for 
what we will call the "Abelian" and for the "non-Abelian" contributions. In particular, 
the structure of the interaction Hamiltonian is independent of the functional dependence 
of / and M and determines the general form of and the anisotropy coefficients appearing 
in the final "non-Abelian" expressions (see Sec. 9). When it comes to explicitate the 
wavefunctions, a choice that can help keeping the result as easy to generalize as possible 
is the following 



6B' 



for the transverse mode and 
55 1 



J- 



3/2 



X) 



3/2 



» 



(8.14) 



(8.15) 

-kr]) [97, 98]. Let us see 



n[x)m^ , 

for the longitudinal mode [n is a unknown function of x \ 
why. As previously stated, for f c:^ a°' and with a = 0, 1, —2, it is possible to verify that 
the (physical) transverse mode behaves exactly like a light scalar field in a de Sitter 
background (see also Appendix 15.4). Considering the solution (8.14) then takes into 
account at least these special cases. As to the longitudinal mode, a parametrization was 
adopted as in (8.15) in order to keep the analysis more general and given that, because of 
the instability issues, introducing this degree of freedom into the theory requires special 
attention. We are going to keep the longitudinal mode "alive" in the calculations we 
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present, by considering a nonzero function n{x), and focus on the simplest case of 
/ = 1. This case is known to be affected by quantum instabihties in the longitudinal 
mode, anyway we choose / = 1 for the sake of simplicity in our presentation. The 
results can be easily generalized to gauge invariant models (please refer to Sec. 12 for a 
sample generalization of some of the calculations to massless / ~ a^^'"^-* models). 
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9. Correlation functions of ( in the SU{2) model 



We are now ready to compute the power spectrum, bispectrum and trispectrum for the 
curvature fluctuations C generated during inflation 

(C^^C,^; = {2ixf5^'\h + k2)Pdk), (9.1) 
{Ck,C>:,M) = i27r)'S^'\h + h + h)Bdku h) (9.2) 

iCkSk^CkM) = (2^)''^^'^(^i + k2 + h + h)Tc{h, h, h, h). (9.3) 

Notice that, on the right-hand side of (9.1) through (9.3), we indicated a dependence 
from the direction of the wavevectors; in models of inflation where isotropy is preserved, 
the power spectrum and the bispectrum only depend on the moduli of the wave vectors. 
This will not be the case for the SU{2) model. 

The (5N formula (3.2) will be applied to our inflaton+S'f/(2)vector model 
C(x , t) = N^5<P + mjAl + \n^^ {5<pf + ^iVjM^M^ + N^dcj^SAl 

+ ^N^H^m' + ^K^c]'5Al5Al5Al5A'i^ + (9.4) 
where now 

and so on for higher order derivatives. 

Our plan is to show the derivation the correlation functions of C, from the ones of 50 
and 5A^, after a replacement of the (5N expansion (9.4) in Eqs. (9.1) through (9.3). 
The correlation functions can be evaluated using the Schwinger-Keldysh formula (2.1), 
that we recall here for convenience 



(i]|e(t)|fi) = (0 



AjlHi{t')dt' 



0). (9.6) 



When calculating the spectra of (, the perturbative expansions in Eq. (9.4) and (9.6) will 
be carried out to only include tree-level contributions, neglecting higher order "loop" 
terms, either classical, i.e. from the 6N series, or of quantum origin, i.e. from the 
Schwinger-Keldysh series. Assuming that the SU (2) coupling Qc is "small" and that we 
are dealing with "small" fluctuations in the flelds and given the fact that a slow-roll 
regime is being assumed, it turns out that it is indeed safe for the two expansions to be 
truncated at tree-level. 



The correlation functions of ( will then result as the sum of scalar, vector and (scalar 
and vector) mixed contributions. As to the vector part, this will be made up of terms 
that are merely generated by the 6N expansion, i.e. they only include the zeroth order 
of the in-in formula (we call these terms "Abelian", being them retrievable in the U{1) 
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case), and by ( "non-Abelian" ) terms arising from the Schwinger-Keldysh operator ex- 
pansion beyond zeroth order, i.e. from the gauge fields self-interactions. 



Let us now discuss the level of generality of the results we will present in the next 
sections, w.r.t. the choice of a specific Lagrangian. 

The expression for the Abelian contributions provided in Sees. 9.1 and 9.2 apply to any 
SU{2) model of gauge interactions with no direct coupling between scalar and vector 
fields (extra terms would be otherwise needed in Eqs. (9.18) and (9.19)). The next 
stage in the Abelian contributions computation would be to explicitate the derivatives 
of the e-foldings number and the wavefunctions of the fields: they both depend on the 
equations of motion of the system, therefore the fixing of a specific model is required at 
this point. 

As to the non-Abelian contributions, the results in Eqs. (9.33) and (9.34) are completely 
general except for assuming, again, that no direct vector-scalar field coupling exists. The 
structure of Eqs. (9.49) and (9.52) is instead due to the choice of a non-Abelian gauge 
group. The expressions of the anisotropy coefficients /„ and L„ in Eqs. (9.49) and (9.52) 
depend on the specific non-Abelian gauge group (for SU{2) one of the In is given in 
Eq. (9.51)). Finally, the specific expressions of the isotropic functions F„ (a sample of 
which is shown in Eq. (9.50)) and G„ were derived considering the Lagrangian (8.11) 
with / = 1 and the eigenf unctions for the vector bosons provided in Eqs. (8.14) and 
(8.15). 



9.1. The power spectrum 

The power spectrum of C, can be straightforwardly derived at tree-level, using 
expansion (9.4), from the infiaton and the vector fields power spectra 

P^ik) = P^'^k) [l + g''^ [k ■ Na) {k ■ N,) + is'^^'k ■ [Na x Nb 
The isotropic part of the previous expression has been factorized in 

P'^\k) = NlP^{k) + (iV, ■ iV,) , 
where we have defined the following combinations 



P. 



ab 



;i/2)(pj 



ab 



R ±Pl 



ab\ 



from the power spectra for the right, left and longitudinal polarization modes 

P^' ^5ab6AUk,t*)6A'*{k,f), 



Pf =6abSAUk,t*)6Ai*{k,t* 



Plong = ^ab5A1^^g{k,t*)5A\l^g{k,t*). 

The anisotropic parts are weighted by the coefficients 

N^N' [Ptg - Pf 



..ab 



NlP^ + (A^, ■ N,) Pf 

jya jyfc pab 

NlP^+(N,-N,)pf 



the (5N 
(9.7) 

(9.8) 

(9.9) 

(9.10) 
(9.11) 
(9.12) 

(9.13) 
(9.14) 
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(where a sum is intended over indices c and d but not over a and b). Eq. (9.8) can also 
be written as 

pcd- 



cd~ 



after introducing the parameter 



(9.15) 



(9.16) 



Notice that what when we say "isotropic" , as far as the expression for the power spec- 
trum is concerned, we simply mean "independent" of the direction of the wave vector. 
In this case instead, the vector bosons introduce three preferred spatial directions: the 
r.h.s. of Eq. (9.7) depends on their orientation w.r.t. the wave vector. 

As expected, the coefficients g"''^ and s"^ that weight the anisotropic part of the power 
spectrum are related to /3cd, i.e. to the parameters that quantify how much the expansion 
of the universe is affected by the vector bosons compared to the scalar field. 
Assuming no parity violation in the model, we have s'^^ = 0; the parameters g"''^ and (5ab 
are instead unconstrained. In the U{1) case and for parity conserving theories, Eq. (9.7) 
reduces to [92] 



P^{k) = P^'^k) [l + g{k 



n 



(9.17) 



where n indicates the preferred spatial direction; also one can check that in this simple 
case, if P+ ~ P^ and Piong = kP+ {k ^ 1), the relation g = {k — 1)13 /{I + (3) holds, 
where (3 = {Na/N^Y (the anisotropy coefficient g is not to be confused with the SU{2) 
coupling constant gc). If it is safe to assume l^?! ^ 1 (see discussion following Eq. (1.1) 
and references [70, 71]), a similar upper bound can also be placed on (3. 
In the case where more than one special directions exists, as in the SU{2) model, no 
such analysis on the anisotropy data has been so far carried out, the g"' parameters 
cannot then be constrained, unless assuming that the three directions converge into a 
single one; in that case a constraint could be placed on the sum 1(71 = 1 J2a9°'\j where 



a = 1,2,3 and P^(A;) 



TDiso 



{k) 



9" 



k ■ hr 



In the next sections we will present the results for the tree-level contributions to the 
bispectrum and to the trispectrum of C,. 

These can be classified in two cathegories, that we indicate as "Abelian" and "non- 
Abelian". The former are intended as terms that merely arise from the 5N expansion 
and are thus retrievable in the Abelian case; the latter are derived from the linear and 
quadratic expansions (in terms of the gauge bosons interaction Hamiltonian) of the 
Schwinger-Keldysh formula and are therefore peculiar to the non-Abelian case. 
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9.2. Bispectrum and trispectrum: Abelian contributions 



By plugging the 5N expansion (9.4) in Eqs. (9.2) and (9.3), we have 
B^ih, h, h) 3 ^N^N^^ [P^{ki)P^{k2) + perms.] 



perms. 



+ -N^KN;, [P^{h)Ilf,{k2)+ perms. 

+ N^NlB^{ki,k2,k3), 
for the bispectrum and 
T^(A;i, k2, ks, k^) D N^T^{ki, k2, ^s, h) 

+ N^N^^, [P^{ki)B^{\ki + k2\,k3, ki)+ perms. 

+ NINI^N;, [p;l{h)B^{k,, k2, \h + h\)+ perms. 

+ NlNl^ [P^{k^)P^{k2)P^{\h + h\)+ perms'. 

+ NlN^^^ [P^{k,)P^{k2)P^{k3)+ perms.] 

+ NIN^N;, [P;l{k, + h)P^{k{)P^{k2)+ perms'. 

+ NI^N^.NI^N;, [p;,;{h)Ptt{k2)P^{\k + k,\)+ perms. 

+ NI^HN;^, [P^{k{)P^{k2)P;l{h)+ perms: 

+ N^Ni:^N^,Ni:, [p;;{k{)Ptt{k2)P^{k,)+ perms. 

+ N^^N^^^N'i; [P^{k2)P^{\h + ^2|)P;'(^4) + perms." 

+ NI^^NPN;,N^ [i^r(fc2)n7(^i + k2)P^{k,)+ perms. 

+ Nj^N^N^^N',] [p;;{k)Ptt{k2)Pf,{h + h)+ perms. 



def 



(9.18) 



(9.19) 



P^'^{k,)P',l{k2)P;:f{k,)+ perms. , 
for the trispectrum. 

Before we proceed with exphcitating these quantities and for the rest of the thesis, the 
A'^ coefficients will be set to zero: it is possible to verify that the temporal mode B^ = 
is a solution to the equations of motion for the vector bosons, after slightly restricting 
the parameter space of the theory (see Appendix 15.4); the adoption of this kind of 
solutions, which is related to the assumption of a slow-roll regime for the vector fields, 
implies that the derivatives of w.r.t. the temporal mode can be set to zero. 
Let us now provide some definition for the quantities introduced in (9.18)-(9.19): we 
are going to switch from the greek indices ^, u, ... to the latin ones, generally used for 
labelling the three spatial directions, in order to stress that all of the vector quantities 
will be from now on three-dimensional 

ng'(fc) ^ T^p-{k)Pl'> + zT^f\k)P^'' + Tlj^\k)P^^^, (9.20) 

where 

\k)^e'^{k)ef{k) + ef{k)ef{k), (9.21) 
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rpeven / 
ij 



e^{k)ef{k)-ef{k)ef{k) 



R( 



(9.22) 
(9.23) 



Tl]-\k)^e\{k)ef{k). 

The polarization vectors are e^{k) = ^(cos^ cos0 — i sin 0, cos ^ sin (p + i cos 0, — sin 9), 
e^{k) = e*^{k) and e\k) = k = (sin ^ cos 0, sin ^ sin 0, cos ^), from which we have 

\k) 



rpeven i 
ij 



hi hj , 



ilong 



{k) 



(9.24) 
(9.25) 
(9.26) 



The purely scalar terms in Eqs. (9.18)-(9.19) are already known from the literature |. 

As to the mixed (scalar-vector) terms, they can be ignored if one considers a Lagrangian 
where there is no direct coupling between the infiaton and the gauge bosons but the 
latter condition is not sufficient for concluding that the mixed derivatives are null. As 
an example, it is useful to refer to [148] which, among other things, includes an analytic 
study for the case of a set of slowly rolling fields with a separable quadratic potential. 
The number of e-folding is written as a sum of integrals over the different fields, to be 
evaluated between their values at an initial (generally set at around horizon crossing) 
and a final times. For each field, the value at the final time depends on the total field 
configuration at the initial time, so the mixed derivative of can in principle be non- 
zero. Anyway, if the final time approaches the end of infiation, it is reasonable to assume 
that, by then, the fields have stabilized to their equilibrium value and no longer carry 
the memory of their evolution. If this happens, the sum of integrals which defines A^ be- 
comes independent of the final field configuration and its mixed derivatives can therefore 
be shown to be zero. It turns out that we are allowed the same kind of analytic study, 
if we work with the Lagrangian in Eq. (8.11) and introduce some slow-roll assumptions 
for the vector fields (see Appendix 15.4 for a discussion about these assumptions and 
Appendix 15.5 for the actual calculation of A^ and its derivatives). 



Let us then look at the (purely) vector part. Its anisotropy features can be stressed by 
rewriting them as follows 



B^ki, k,, ks) D -ArXAr^^jn-(/ci)n5f (/C2) = M.^AT^^jMf 

T^{k„ k„ h, h) D N^N;:N^^N'^]p^;^{h)P'J{h)Pf,iki + h) 

+ A^.^A^r We';P;.'(^l)P.1(fc2)P,t(4) 

= M^LgM^ + Mi M^MtNf,, 

where 

Ml{k) = NlPt,%k) = Pl\k) [5,uNl+p''%k)kk {k ■ N,) + zq^^k) {k x 



(9.27) 



(9.28) 



(9.29) 



X In single-field slow-roll inflation = H^/2k^, where is the Hubble rate evaluated at horizon exit; 
the bispectrum and the trispectrum of the scalar field (B^ and T^) can be found in [12, 13, 147, 152. 39]. 
For the bispectrum see also Eq. (3.10) 
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In the previous equations, we defined 

pac 



q^\k) 



pac _ 
long 

Pi 

pac 



tq'f{k)k-Ni,xNl^]. (9.30) 

(9.31) 
(9.32) 



with = {Nl Nl Nl) and Ni, ^ {NH iVf , nH). 

Notice that, as for the power spectrum (9.7), also in Eqs. (9.27)-(9.28) the anisotropic 
parts of the expressions are weighted by coefficients that are proportional either to 
P_ or to {Piong ~ P+)- When these two quantities are equal to zero, the (Abelian) 
bispectrum and trispectrum are therefore isotropized. P_ = in parity conserving 
theories, hke the ones we have been describing. According to the parametrization (8.15) 
of the longitudinal mode, we have Piong — P+ = {\n{x)\'^ — 1)P+. 

9.3. Bispectrum and trispectrum: Non-Abelian contributions 
We list the non-Abelian terms for the bispectrum 

B^ku h) 3 KNlN^Bt^^ih, h) (9.33) 
and for the trispectrum 

T^ih, h, h, h) D KNiN^N',Tff,f{h, h, h, h) 

+ NiNiN^Nl [p^(fc3)SS(fci, ^2, h + + perms'. 

+ KNIn'^N'J: [P^\k,)Bf^^{k, + ^2, 4, k,) + perms. 

The computation of the vector bosons spectra 

{5At5A]5Al) = 5(3) (^^ + ^ + k,)B^^^, 

{6A16A'^6Al6Af) = 6^'\h + h + h + h)T$f, 

will be reviewed in this section. This requires the expansion of the in-in formula up to 
second order in the interaction Hamiltonian 



(9.34) 

(9.35) 
(9.36) 



(9.37) 

+ ^(T[e 1^ dr^' [HUri) - HUr]')) J'^ drj' [HtM) " HUrl')) ])■ 

The interaction Hamiltonian needs to be expanded up to fourth order in the fields 
fluctuations, i.e. Hint = H-^l + Hinh where 

i/g) _ g^£-'>^g'''g^^ [diSB]) SB16B^ + e'"^ q'^ B"} 5 B]5 B15 Bf (9.38) 

h\S, = gy^'e'^'g^^g'^dBtSBlSB^^dBf. (9.39) 

To tree-level, the relevant diagrams are pictured in Figs. 5 and 6. By looking at 
Eqs. (9.38) and (9.39), we can see that there is a bispectrum diagram that is lower 
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Figure 5. Diagrammatic representations of the tree- level contributions to the vector 
fields bispectrum. 




Figure 6. Diagrammatic representations of the tree- level contributions to the vector 
fields trispectrum: vector-exchange (on the left) and contact-interaction (on the right) 
diagrams. 

in terms of power of the SU{2) coupling (~ Qc) compared to the trispectrum (~ g1)\ as 
a matter of fact, for symmetry reasons that we are going to discuss later in this section, 
interaction terms are needed to provide a non-zero contributions to the bispectrum. 



The propagators for "plus" and "minus" fields are 

5Bt^+{ri)5B]^^{r^") = mv\v"Mv' " v") + n^^iv'.vMv" " V), (9-40) 

5i?r(r/>4-(r/") = n^^(V, V), (9.41) 

6Br^B^'^^") = n^nV, V), (9.42) 

SBrivlSB'^'^iv") = n^;{v\v")e{v' - v") + t^tHv'^v'Mv" - v'), (9.43) 

or 

ng'(fc) = 7;7"(fc)Pf + iT°f\k)P^I> + Tlp'{k)P,f (9.44) 

ng'(fc) = 7;7^"(^)Pf + iT°f^{k)P^I> + Tlp'{k)P,f (9.45) 



in Fourier space. In the previous equations we set P^^ = {1/2){P^ ± P^'^), P^ = 
5ab5B'^{k,r]*)SB*£\k,ri) and Pf = (Pf )* (similar definitions apply for Pf and P;f„^). 

We are now ready to show the computation of the following contributions to the 
bispectrum and trispectrum of ( 

{C>^Sk.C>^,) ^ KNlN^{5Ank\)SA%k2)SAl{ks)), (9.46) 

iQSk,Ck,CkJ ^ KNlN^N!,{6A^{h)SA%h)6Al{ks)6Af{h)). (9.47) 
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Eq. (9.46) becomes 



kS^^Kh + k2 + h] 



perms. + c.c. 



(9.48) 



Even before performing the time integration, one realizes that, because of the 
antisymmetric properties of the Levi-Civita tensor, the ~ Qc contribution on the r.h.s. 
of Eq. (9.48) is equal to zero once the sum over all the possible permutations has been 
performed. The vector bosons bispectrum is therefore proportional to g^. The final 
result from (9.48) has the following form 



{Ck,Ck,Ck,) ^ i27rfS^^\h + k2 + h)9!H^T.Fn{h,vlUk-k^A,-A,,k-A,) (9.49) 



n 



where the F^s are isotropic functions of time and of the moduli of the wave vectors 
{i = 1,2,3) and the /„'s are anisotropic coefficients. The sum in the previous equation 
is taken over all possible combinations of products of three polarization indices, i.e. 
n e {EEE,EEl,ElE,...,lll), where E stands for "even", I for "longitudinal". The 
complete expressions for the terms appearing in the sum are quite lengthy and we list 
them in Appendix 15.6. As an example, we report here one of these terms 

Em = - ^^(^*) 24fc6fc2^2^2^,.2 ["^EEE + {Beee cosx* + CEEESmX*) EiX*] (9.50) 

Iiu = (A;i ■ N'^') (k ■ N'') {k ■ N'') {k ■ k) {k ■ A') 

- [k ■ N'' ) (k2 ■ N'') (ki ■ N'') [h ■ k) [k ■ i^) ) + (1 ^ 3) + (2 ^ 3)] (9.51) 

where Aeee, Beee and Ceee are functions of x* and of the momenta ki = \ki\ (they 
are all reported in Appendix 15.6), Ei is the exponential-integral function. As we will 
discuss in more details in Sec. 11.2, one of the more interesting features of these models 
is that the bispectrum and the trispectrum turn out to have an amplitude that is mod- 
ulated by the preferred directions that break statistical isotropy. 

Let us now move to the trispectrum. Again, we count two different kinds of 
contributions, the first from ~ g^. and the second from ~ g"^ interaction terms, 
respectively in iff^] and H^^^. The former produce vector-exchange diagrams, the latter 
are represented by contact-interaction diagrams (see Fig. 6). Their analytic expressions 
are different, but they both have a structure similar to (9.49) 

(CkSkMCk,) ^ i27r?6^'\h + h + h + h)glHl (9.52) 
X Gn{ki, ki2i ^14) ^7 )Ln{ki ■ kj, Ai ■ Aj., ki ■ Aj) 

n 

where we define ky^ = \ki + k2\ and /cf4 = j/ci + A;4|. We will present the details of the 
computation of the gauge fields trispectrum and the explicit expression of the functions 
appearing in (9.52) in the next section. 



42 



9.3.1. Trispectrum from vector bosons: exchange diagram 

Let us begin with the two vertex diagram. Using the language of Eq. (9.37), it can be 
put in the form 

(e(r/*)) D ^(T[e f_^Jri {h+ {r^ ) - {r^ )) J'jv" {h^v") - H~ (v"))]) (9.53) 

where now H = H- '^il , 6 = and the inclusion symbol as usual points out 

that what stands on the right-hand side is only one of the contributions to (6(77*)). 
Eq. (9.53) can be rewritten as follows 

(e(r/*)) D ^^(T [9(^ + 5+ C+D)]) (9.54) 



where 



A = r drlH+ijI) r dri"H+{r]") 
J— 00 J— 00 

B= r driH~{ri) T d7]"H-{r]") 
J— 00 J— 00 

C=- r d7]'H+{r]') r dr]"H~{7]") 

D=- r d7]'H^{r]') r dr]"H+{7]"] 

J —00 J— 00 



For each one of the integrals listed above, due to the presence of both the fields and their 
spatial derivatives in H^J^, there are three different sets of contractions of the external 
with the vertex field-operators: for the first set, the field-operators with derivatives in 
the vertices are both contracted with external fields; for the second one, only one of the 
two field-operators with derivatives contracts with an external field (the other contracts 
with another internal field); for the third set, the field-operators with derivatives contract 
with each other. 

A sample set of contractions of the first type is provided in the following equation 

!, J 0^ '1.' ' "i," " 

rpabcd He o c a o c r, 

J dr^'a^r]') J dr^' a\7i')g^^g^ ^' g-'^g^ ^j' (9-55) 

where the first four lis correspond to contractions between external and internal fields 
whereas the last one indicates the contraction between the two remaining internal 
fields. The factor comes from expressing the external (physical) fields in terms 
of the comoving ones. As a reminder, we define (27r)^5*^^^(fci + k2 + + k4)T^I'^f = 
{5A'^5A'^j5A15A'l). The expression in (9.55) can be rewritten as follows 



X 



I t t II tt ft I I It It I It 

•abed JC ^a b c ^a b e 



rpaoea 



_^at,e^a b e hh,§aa ^bb ^ea ^dh ^e e 



X 



(/ / 



afi-yScr 



T^OLrjijB rT-17 rjiS rjKT 
m jm im nn 



(9.56) 
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where the greek indices of the sum indicate either the transverse (E) or longitudinal (1) 
modes, the (^J drj' J drj"^ stand for the integrals over the wave functions, the chosen time 

variable being x = —kr] {k = J2i=i,...,4ki, ki = \ki\). 

Let us define the coefficients T°:l7^" = kmk'T'^T^^T?',Tf,T'',. They should be 
calculated for each one of the three different sets of contractions and for each permutation 
within the specific set. This is a straghtforward but rather lengthy and not particularly 
interesting calculation. A convenient way to proceed could be the following: we first 
compute the time integrals in order to find out which one among the combinations 
of longitudinal and transverse mode functions in the string [a, /3, 7, 5, a] provides the 
highest amplitude for the trispectrum (in order to be able to perform this comparison we 
work, as it is usually done when trying to quantify the amplitude of a three or of a four- 
point function, in the so called "equilateral configuration", which for the trispectrum 
means taking k\ = k^ = kj, = k^\ for the combination with the highest amplitude, we 
then calculate the coefficients T°!j^^^" for all the different sets of contractions and sum 
over all the permutations. 

Let us now perform our calculations. The wavefunctions we are going to adopt were 
introduced in Sec. 8.5 (see Eqs. (8.14) and (8.15)). It is possible to verify that B = A* 
and D = CT and that integrals of type A are consistently smaller in amplitude than 
integrals of type C. We therefore report the combined contribution C + D = 2Re[C\ for 
one of the permutations 



dt] / d7] 



EEEEE 



8klk^kiklkl~^{k-f2 + ki + k2){k^2 + h + kA)x 
X [(M - 2E)[{N - 2F){AB + CD) + {2H + L){CB - AD)] 
+ {2G + P)[{N - 2F){AD - CB) + (2H + L){AB + CD)]] 



(9.57) 



df] 
dr] 
dr] 
dr] 
dr] 



dr] 
dr] 
dr] 
dr] 
dr] 



EEEEl 



EEEll 



EElll 



EUU 



UUl 



n [x 



4 / * 

n [x 



n^ix* 



n^ix* 



dr] 
dr] 
dr] 
dr] 
dr] 



dr] 
dr] 
dr] 
dr] 
dr] 



EEEEE 



EEEEE 



EEEEE 



EEEEE 



EEEEE 



(9.58) 
(9.59) 

(9.60) 
(9.61) 
(9.62) 



where A, B, C, D, E, F, G, H, L, M, N and P are functions of x* and of the momenta 
moduli to be provided in Appendix 15.7 (see Eqs. (15.7.1) through (15.7.9)). Obviously, 
the value of the integrals does not change when permuting its labels aP'~f6a, apart from 
a different power of the coefficient n{x*). We need now to find out if there is one, 
among the integrals in Eqs. (9.57) through (9.62), that has the largest amplitude, i.e. 
understand if something can be said about the order of magnitude of n{x*). We could 
try to extrapolate some information about n{x*) from what happens at very late times. 
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In the models discussed in Ref. [135, 92] it turns out that the longitudinal mode is 
5-8" = \/25B^ §. If this is the correct asymptotic behaviour and we find it reasonable 
to extrapolate back until the horizon crossing epoch, it is then correct to conclude 
that in this case the amplitude is the largest for the integral among the ones listed in 
Eqs. (9.57) through (9.62) containing the highest powers of n, i.e. for (^J drj' J dr]"^ 
The coefficients we intend to calculate are then of the kind Tj^jH only. We list them below 
for the three different sets of contractions we introduced above and for one particular 
permutation (see again Appendix 15.7 for more details) 



rpiiim) ^ kikskusA {h ■ k^2) {ki ■ h) [h ■ h) [h ■ ^12) , (9-63) 
rpiiim ^ k3kf2ki234 {h ■ ^12) (^2 ■ ^12) (^3 ■ ^4) [h ■ ^12) > (9-64) 
rpims) ^ k{^k^^ki2u {h ■ ^12) (^2 ■ ^12) (^3 ■ ^12) (^4 ■ ^2) • (9-65) 

We adopted the following notation: kg = \ks\, kg = kg /kg, the index s running over 
the four external momenta; kss's"s"' = ksiks'jksiikksini^ with s,s',s", s'" = 1,2,3,4 and 
with the indices j, fc, / indicating the spatial components of the vectors; k^y = ks + kg', 
k^y = \ks + ks'\ and so k^^, = ky^,/ky^,. 

It is possible to prove that, once the Levi-Civita coefficients and the sum over the 
permutations are taken into account, only the first set of contractions is left. The final 
result after these cancellations can be written in the following form 



{SA^SA'^SAlSAf), D {27if6^^\h + k2 + h + h)gl e'"^' [l x A;i234 x (jy, 

+ IIX fci324 X (j2 + III X ^1432 X (j^ ] ' (9-66) 

All the possible permutations have been included in the previous equation and, as a 
reminder, the indices i,j, k, I are hidden in kss's"s"' on the right-hand side. We define 

/ = X ^ ^ ^ 



^8k'fk^k'iklkl-^{kf2 + ki + k2){kf^ + ^3 + ki)x*^ J 
X [(M - 2E) [{N - 2F) {AB + CD) + {2H + L) {CB - AD)] 

+ {2G + P)[{N - 2F){AD - CB) + {2H + L){AB + CD)]] (9.67) 

(from Eqs. (9.57)-(9.62)). The function // is defined from / by exchanging /c2 with k^ 
and kf2 with k-f^; similarly, /// is defined from / by exchanging k2 with ^4 and kf2 with 
/cf4, so they are all functions of the horizon crossing time x* = —krj* and of the moduli of 

the external momenta and of their sums. This amounts to seven independent variables, 

pa/37. 

■ ijkl 



X*, ki, /c2, /c3, ki, kf2 and k-f^. The coefficients U {i = 1,...,12) come from T^^i^'^ 



§ As another example, in models with varying kinetic function and mass [95, 96], we have verified that 
n{x) ^ 1 at late times and for a vector field that is light until the end of infiation. 
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so they are also functions of the momenta moduh (see Eqs. (15.7.10) through (15.7.21) 
for their expressions). Finally, the anisotropic part of Eq. (9.66) is represented by the 
kss's"s"' terms, which, in the final expression for the curvature perturbation trispectrum, 
have their spatial indices contracted with the derivatives N°- of the number of e-foldings 
w.r.t. the vector fields as follows 



{Q^Ck^CksCkJ ^ N^N^mfmSA'^SAlSAf). D (2vr)35(3)(fci + h + h + k,)gl 

X J X tj^ X A/ + perms. 

where the anisotropic term in the first permutation is 

A/ = e"^''" e""^"" N^N!lN^Nfh234. 



(9.68) 



(9.69) 



It can be interesting to rewrite A/ in terms of all its variables 



{a<b)a,b=l 



(9.70) 



(iV")2(^iV^) X n f^et (m)'^'"'') 

[ij] = [l,2],[3,4] 

The M/'s are 2x2 matrices whose entries are represented by the cosines of the angles 
between the wavevectors and the N"- 



COS 6ia COS 6ja 

COS 9ib COS 9jb 



i.e. cos6'ja = ki ■ N"' and so on. 

The two permutations in Eq. (9.68) can be written in a similar fashion with anisotropic 
coefficients 



[ij] = [l,3],[2,4] 
[i,i] = [l,4],[3,2] 

The number of angular variables is equal to 12. These are to be added to the six 
scalar variables from the isotropic part of (9.68) (/ci, k2, k^, k^, k^2 ^"^^ ^14) to 
three parameters represented by the lengths of the vectors A^'* in A. The anisotropy 
coefficients here become equal to zero in the event of an alignment of the gauge vectors 
along a unique direction. 



Ajj = e''^' e"^ N;'N^N^Nfkn24 = E 

{a<b)a,b=l 



Ajjj = e"^^ e^'^^ N^N!;N'^Nfku32 = E 

{a<b)a,b=l 



9.3.2. Trispectrum from vector bosons: point-interaction diagram 
Let us now move to the one-vertex diagrams 







dr^'(H^r^')-H-{r^')) ) 



(9.71) 
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where now H = if^^] and again 6 = 5A'^^5A\,5A'p5A'^^. After working out the Wick 
contractions, this becomes 



Y.(fdx) T^^T^^TIT'^^ + permutations. (9.72) 

-..R^.x V-' / 0/3-7(5 



X 

Let us hst the coefficients T^^l^ = T^^T^jT^'^Tpj for one of the permutations 



rpEEEE 
mnop 


^mo^np ^mokpikfiA 


- 5mokn2kp2 + 5mokn2kpik2 " /i;4 ^npkozkmZ 






+ /C3434 + ^3232 — ^2 ■ ^4 (^3234 + ^1214) — ^npkmlkol + ki4i4 + /i;i212 






+ Snpkmikoski ■ k% ki 


■ h {ki2-i2 + A;i434) + A;i ■ A;3A;2 ■ A;4, 


(9.73) 


' 1 1X1' il il t 

mnop 


— ^mokpAkni — Smokn2kp4,k2 " ^4 ~ ^3434 + k2 ■ k^ (/i:3234 + ^1214) ~ ^1414 






+ kiis^ki ■ ks — ki2S4ki 


■ A;3A;2 ■ ki, 


(9.74) 


/ 1 T 

mnop 


— SnpkoskmS — ^npkmlkoski ■ k^ — /C3434 + ki ■ k^ (/i;i434 + /C1232) — ^3232 






+ ^3234^2 ■ ^4 ~ ^1234^1 


■ ■ ki, 


(9.75) 


rjiElE E 

mnop 


— Smokn2kp2 — Smokn2kp4,k2 ■ ki — /C3232 + /^2 " ^4 (^3234 + ^1214) — ^1212 






+ ^1232^1 ■ ^3 ~ ^1234^1 


■ hk2 ■ ki, 


\\). Co j 


rjilE E E 
mnop 


— kmlkolSnp — kmlkoskisSnp — /Cl414 + ^1 " ^3 (^1434 + ^1232) — ^1212 






+ ^1214^2 ■ ^4 — ^1234^1 


■ hh ■ ki, 


f9 77) 


rpEEU 

mnop 


= /C3434 — ^3234^2 ■ ^4 — 


kiisiki ■ ks + /ci234^i ■ ^3^2 ■ ki, 


(9.78) 


rpElEl 

mnop 


= 5mo^n.2^p4^2 ■ ki — k2 


■ ki (A;3234 + ^1214) + ^1234^1 " ^3^2 " ki, 


(9.79) 


rpEUE 

mnop 


= ^3232 — ^3234^2 ' ^4 " 


^1232^1 ■ h + ^1234^1 ■ A;3A;2 ■ ki. 


(9.80) 


rpUEE 

mnop 


= ^1212 — ^1214^2 ■ ^4 — 


ki2?.2ki ■ h + A;i234A;i ■ A;3A;2 ■ ki. 


(9.81) 


rplEEl 

mnop 


= A;i4i4 — /Cl214^2 ■ ^4 — 


kiisiki ■ ks + /i;i234^i ■ ^3^2 ■ ki. 


(9.82) 


rplElE 

mnop 


= Snpkmikoski ■ k^ — ki 


■ k^ (A;i434 + ki2^2) + ki23iki ■ k2,k2 ■ ki. 


(9.83) 


rplllE 

mnop 


= ^1232^1 ■ ^3 ~ ^1234^1 


■ hh ■ ki. 


(9.84) 


rpUEl 

mnop 


= ki2iik2 ■ ki — A;i234^i 


■ hk2 ■ ki. 


(9.85) 


rplEU 

mnop 


= kiisiki ■ ks — A;i234^i 


■ hh ■ ki. 


(9.86) 


rpElU 

mnop 


= ^3234^2 ■ ^4 — ^1234^1 


■ hh ■ ki. 


(9.87) 


rpllll 

mnop 


= ^1234^1 ■ ^3^2 ■ ^4- 




(9.88) 


When evaluating the integrals (/ dx)^^^^, we use again the wavefunctions 


previously 



introduced in Eqs. (8.14) and (8.15). The final result is 



dx] = 0,2,2/2 *7 [Qeeee + Aeeee cix* {Beeee cos X* + Ceeee sin x*) 

J EEEE 24:k^k(k2k^kfx*' 

+ Deeee six* {Eeeee cosx* + Feeee sinx*) ], (9.89) 
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dx = n\x*) [ dx] , (9.90) 

/ EEEl \J / EEEE 

dx) = n\x*) ( f dx) , (9.91) 

/ EEll \J / EEEE 

dx) = n^{x*) ( f dx) , (9.92) 

/ lllE \J ) EEEE 



dx = n%x*) { dx) , (9.93) 

/ nil \J J EEEE 

where Qeeee, Aeeee, Beeee, Ceeee, Deeee, Eeeee and Feeee are functions of x* 
and of the momenta ki = \ki\, ci and si stand respectively for the Coslntegral and the 
Sinlntegral functions. The expressions of these functions can be found in Appendix 15.8. 
It is again important noticing that the anisotropy coefficients become zero if the gauge 
fields are all aligned. 

Finally, summing up the coefficients in Eqs. (9.73) through (9.88), one realizes that if 
the longitudinal and the transverse mode evolve in the same way, the total contribution 
from the point-interaction diagram is isotropic 

^ ^j IQeEEE + AeeEECIX* {BeeEE COS X* + CEEEE^mX*) 

+ Deeeesix* {Eeeee cosx* + Feeee sinx*) ] + permutations. (9.94) 
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10. Amplitude of non-Gaussianity: /tvl and tnl 

Our definitions for tlie non-Gaussianity amplitudes are 



6 Eg(fci,fc2,fc3) , . 

^JNL pisoi^ki)P'''°{k2)+ perms. ^ 

piso(^ki)P''°{k2)P''°{kf^) + 23 perms. ^ ' ' 

The choice of normalizing the bispectrum and the trispectrum by the isotropic part of 
the power spectrum, instead of using its complete expression P^, is motivated by the fact 
that the latter would only introduce a correction to the previous equations proportional 
to the anisotropy parameter g, which is a small quantity. 

The parameters Jnl and tnl receive contributions both from scalar ("s") and from 
vector ("w") fields 

I'nl = Inl + /nx) (10.3) 
r^L = ra + rS. (10.4) 

The latter can again be distinguished into Abelian (A) and non-Abelian (NA) 

Al = Al' + /K"', (10.5) 

Ai = Al-' + r'„f + A™ ■> + r^"- (10.6) 

The contribution comes from Eq. (9.27), f^i^^ from (9.49), T^l^ and t^^^ from 
(9.28), finally r^^^'^ from (9.52) and r^^^'^ from the last line of (9.34). 

In order to keep the vector contributions manageable and simple in their structure, all 
gauge and vector indices will be purposely neglected in this section and so the angular 
functions appearing in the anisotropy coefficients will be left out of the final amplitude 
results. This is acceptable considering that these functions will in general introduce 
numerical corrections of order one. Nevertheless, it is important to keep in mind that 
the amplitudes also depend on the angular parameters of the theory. 
We will now focus on the dependence of Jnl and tnl from the non- angular parameters 
of the theory and quickly draw a comparison among the different contributions listed in 
Eqs. (10.3) through (10.6). 

The expression of the number of e-foldings depends on the specific model and, in 

particular, on the mechanism of production of the fiuctuations. Two possibilities have 

been described in Sec. 8. For "vector infiation" we have 

A°- ■ S hS^^ 

N' = — ^ N'{ = °^ (10 7) 
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Table 1. Order of magnitude of fpfL in different scenarios. 
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J NL 
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Table 2. Order of magnitude of the vector contributions to tnl 
in different scenarios. 
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v.curvaton 



same as above 
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Table 3. Order of magnitude of the ratios fp^^/ fffi in different scenarios. 
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(see Appendix 15.5 for their derivation). In the vector curvaton model the same 
quantities become [92, 97] 

AT^ = Ir \ , m = ir^^^^. (10.8) 

Neglecting tensor and gauge indices, the expressions above can be simplified as Na — 
A/m?p and N^a — l/?np in vector infiation, Na — r/A and Naa — t/A^ in the vector 
curvaton model. Also we have Naaa = in vector infiation and Naaa — r /A^ m vector 
curvaton. 

We are now ready to provide the final expressions for the amplitudes: in Table 1 we list 
all the contributions to /atl. Table 2 includes the vector contributions to tatl, the scalar 
contributions being given by 

2 

rj;| = ^ + (10.9) 

In the expressions appearing in the tables, numerical coefficients of order one have not 
been reported. Also, m is by definition equal to mp in vector infiation and to A/y/r 
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Table 4. Order of magnitude of the ratios T^i^/rffj^ in different scenarios. 
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in the vector curvaton model; ~ (mpy/e) ^ and A'^^^ ^ '^p^, with e = (0^) / {2rn?pH'^). 

The quantities involved in the amplitude expressions are (3, r, e, Qc, rrip/H, A/rrip 
and A/H. We already know that g and /3 are to be considered smaller than one (see 
discussion after Eq. (9.17)). Similarly, as mentioned after Eq. (8.2), r has to remain 
small at least until inflation ends so as to attain an "almost isotropic" expansion. The 
slow-roll parameter e and the SU{2) coupling are small respectively to allow the infia- 
ton to slowly roll down its potential and for perturbation theory to be valid. The ratio 
irtp/H is of order 10^ (assuming e ~ 10~^). Finally, A/mp and A/H have no stringent 
bounds. A reasonable choice could be to assume that the expectation value of the gauge 
fields is no larger than the Planck mass, i.e. A/mp < 1. As to the A/H ratio, different 
possibilities are allowed, including the one where it is of order one (see Sec. 6 of [97] for 
a discussion on this). 

Let us now compare the different amplitude contributions. The ratios between scalar 
and vector contributions are shown in Table 3 for the bispectrum and Table 4 for the 
trispectrum. We can observe that the dominance of a given contribution w.r.t. another 
one very much depends on the selected region of parameter space. It turns out that it 
is allowed for the vector contributions to be larger than the scalar ones and also for the 
non-Abelian contributions to be larger than the Abelian ones. This is discussed more 
in details in Sec. 6 of [97]. An interesting point is, for instance, the following: ignoring 
tensor and gauge indices, the ratio gcA/H, that appears in many of the Tables entries, 
is a quantity smaller than one; if we consider the different configurations identified by 
gauge and vector indices, we realize that this is not always true, in fact the value of 
this ratio can be ^ 1 in some configurations (see also Appendix 15.4 and Eqs. (15.4.9) 
through (15.4.11) in particular). 

Finally, it is interesting to compare bispectrum and trispectrum amplitudes (see 
Table 5). Again, it is allowed for the ratios appearing in Table 5 to be either large or 
small, depending on the specific location within the parameter space of the theory. For 
instance, the combination of a small bispectrum with a large trispectrum is permitted. 
The latter is an interesting possibility: if the bispectrum was observably small, we 
could still hope the information about non-Gaussianity to be accessible thanks to the 
trispectrum. 

Another interesting feature of this model is that the bispectrum and the trispectrum 
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Table 5. Order of magnitude of the ratios r^^,/ (/jvl ) different scenarios. 
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depend on the same set of quantities. If these correlation functions were independently 
known, that information could then be used to test the theory and place some bounds 
on its parameters. 
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11. Shape of non-Gaussianity and statistical anisotropy features 

Studying the shape of non-Gaussianity means understanding the features of momentum 
dependence of the bispectrum and higher order correlators (see e.g. [150]). If they also 
depend on variables other than momenta, it is important to determine how these other 
variables affect the profiles for any given momentum set-up. This is the far as 

the bispectrum and the trispectrum of the gauge fields are concerned, given the fact 
that they are functions, besides of momenta, also of a large set of angular variables (see 
Eqs. (9.49) and (9.52)). 

11.1. Momentum dependence of the bispectrum and trispectrum 

We show the study of the momentum dependence of the F„ and G„ functions in 
Eqs. (9.49) and (9.52) first and then analyze the angular variables dependence of the 
spectra, once the momenta have been fixed in a given configuration. A natural choice 
would be to consider the configuration where the correlators are maximized. 
The maxima can be easily determined for the bispectrum by plotting the isotropic 
functions F„ and Gn in terms of two of their momenta. These plots are provided in 
Fig. 7, where the variables are X2 = k2/ki and X3 = k^/ki. Each one of the plots 
corresponds to one of the isotropic functions appearing in the sum in Eq. (9.49). It is 
apparent that the maxima are mostly located in the in the so-called local region, i.e. for 
/ci ~ ^2 ^ ^3; three out of the eight graphs do not have their peaks in this configuration 
but, at the same time, they show negligible amplitudes compared to the "local" peaked 
graphs. 

The situation is much more complex for the trispectrum, being the number of momentum 
variables larger than three (fci, ^2, ^3, ^4, ^12 ^ind kf^). The momentum dependence 
of the isotropic functions can be studied by selecting different configurations for the 
tetrahedron made up by the four momentum vectors, in such a way as to narrow 
the number of independent momentum variables down to two. A list of possible 
configurations was presented in [44]. We consider two of them, the "equilateral" and 
the "specialized planar" . 

In the equilateral configuration the four sides of the tetrahedron have the same length 
(fci = /c2 = /c3 = /C4), therefore x = k^2/ki and y = k^^/ki can be chosen as variables for 
the plots. The plots of the isotropic functions of contact interaction and vector exchange 
contributions are provided in Fig. 8. The former (c.z.) shows a constant behaviour in 
this configuration, being independent of /cf2 k^~^. The latter (f.e.(/), v.e.{II) and 
v.e.{III)) diverge as k^^^ {i = 1,2,3 respectively for the three plots) in the limit of a 
fiat tetrahedron, i.e. (k-^-fki) — )■ 0. 

In the specialized planar configuration, the tetrahedron is fiattened and, in addition to 
that, three of the six momentum variables are set equal to one another [ki = k^ = kf^); 
this leaves two independent variables, which can be x = /c2/^i and y = k^/ki. There is 
a double degeneracy in this configuration, due to the fact that the quadrangle can have 
internal angles larger than or smaller /equal to tt, as we can see from the plus and minus 
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1.0 -rm 




Figure 7. Plot of r„ = Q{x2 — X3)Q{x3 — 1 + X2)x|a;|i?„(a;2, X3), where 
we define i?„ = k\Fn. The Heaviside step functions Q help restricting 
the plot domain to the region (0:2, ^3) that is allowed for the triangle 
fci + ^2 + fca = (in particular, we set < 2:2). We also set a;* = 1. 
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v.e.(I) 




Figure 8. Plots of the isotropic functions appearing in the vector fields trispectrum 
(from Eq. (9.52)): c.i. is the contribution from contact-interaction diagrams, v.e.{I), 
v.e.{II) and v.e.{III) are the contributions from the vector-exchange diagrams. The 
equilateral configuration has been considered in this figure. 



signs in the expressions for k-f2 and k^^ [44] 

^ = \/l + ^±f V^-a:^)(4-l/^), (11.1) 

^ = \Jx' + y'-^Tf^i^-x-){A-y-l (11.2) 

The two cases are plotted in Figs. 9 and 10. Notice that divergences generally occur as 
X, ?/ — )■ 0, as X -7- ?/ and (x, y) — t- (2, 2). 

11.2. Features and level of anisotropy 

Statistical homogeneity and isotropy are considered characterizing features of the CMB 
fluctuations distribution, if one ignores the issues raised by the "anomalous" detections 
we presented in the introduction. 

Homogeneity of the correlation functions equates translational invariance and hence 
total momentum conservation, as enforced by the delta functions appearing on the left- 
hand sides of Eqs. (9.1) through (9.3). This invariance property can then be pictured as 
the three momentum vectors forming a closed triangle for the bispectrum and the four 
momenta arranged in a tetrahedron for the trispectrum (see Fig. 11). 
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Figure 9. Plots of the contact interaction and of the vector-exchange contributions 
in the specialized planar configuration (plus sign). 



Statistical isotropy corresponds to invariance w.r.t. rotations in space of the momentum 
(for the power spectrum) and of the triangle or tetrahedron made up by the momenta, 
respectively for the bispectrum and the trispectrum. This symmetry can be broken, 
as it for example happens in the SU{2) case, by assuming the existence of preferred 
spatial directions in the early universe that might be revealed in the CMB observations. 
When this happens, the correlation functions are expected to be sensitive to the spatial 
orientation of the wave number or of the momenta triangles and tetrahedrons w.r.t. 
these special directions. Analitically, the bispectrum and the trispectrum will depend on 
the angles among the vector bosons and the wave vectors (besides the angles among the 
gauge bosons themselves), as shown in the coefficients /„ and L„ appearing in Eqs. (9.49) 
and (9.52). This implies that both the amplitude and the shape of bispectrum and 
trispectrum will be affected by these mutual spatial orientations. The modulation of 
the shapes by the directions that break statistical anisotropy was discussed with some 
examples both for the bispectrum and the trispectrum in our papers [97, 98]. These 
examples are here reported in Figs. 12 and 13. 

In Fig. 12 we show the plot of the vector contribution to the bispectrum of (, properly 
normalized in the configuration 

N3 = Na{0,0,1) (11.3) 
N^ = N^ = NA{sm9cos(f),sm9sm4),cos9), (11.4) 
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where, the (x, z) coordinate frame is chosen to be ^3 = x and ki = k2 = z and 6 is 
the angle between Ni^2 and ks. The coefficients /„ in this configuration become 
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20 - 24 cos 5 + 2 cos 6-12 cos^ 5 + 12 cos^ 6-2 cos^ ^ + 6 cos ^ cos^ 6 



— 2 cos 6 cos 6 + 2 cos 6 cos 6 



4 cos 5 

4 - 2 cos ^ - 6 cos^ 6- A cos^ 5 

- 4 cos^ 6 

- 4 cos^ 5 

- 2 cos^ 5 

4-4 cos^ ^ - 8 cos^ 5 
2 + cos 6* - 3 cos^ 6 - A cos^ 5 
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(11.6 
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^11.10 
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where = {A)/{Na)i A being the background value of the A^'s evaluated at horizon 
crossing. The analytic expression of the 'non-Abelian' bispectrum normalized to the 
ratio {glH'^w? / {k\x2x\) , as a function of the angles 6 and 5 and for fixed values of 
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Figure 11. Representation of momentum conservation for the bispectrum (the 
three momenta form a closed triangle) and for the trispectrum (the momenta form 
a tetrahedron). 



X*, X2 and xs is 

B(^{e,5) ~ — ^[cos2(5(8cos^- 1.4 x 10^) + 3 cos 5(cos^ ^ - cos^ ^ - 11) 

- 11 COS 25 -40 - 6 cos 2^ - cos^(3cos^^ - 30 cos ^ - 10)1, (11.13) 

where we set x* = 1, while X2 and x^ were chosen in the 'squeezed' region, X2 = 0.9 and 
X3 =0.1. 

In Fig. 13 we provide a similar plot, but for the trispectrum from vector-exchange 
contributions and in a different configuration 

N2-k = 0{i = l,..A) 

Ni-ki = cos 6, Ni-k2 = 

N3-k2 = cose, N3-ki = 0. (11-14) 

In addition to that, let us assume that all the have the same magnitude Na- In this 
configuration, we have 

Ai = Am = N^cos'^ecos'^S, A// = 0, (11.15) 

therefore the the expression in Eq. (9.68) becomes 

D g^^H^N^ [ISO] cos^ 6 cos^ 6 (11.16) 

where the expression in brackets includes an isotropic term (which is rotationally 
invariant) 

ISO^I^^^' (^IpU + IIlf:t,y (11.17) 
Fig. 13 plots the trispectrum contribution in Eq. (11.16) normalized to its isotropic part. 

Another comment should be added concerning statistical anisotropy in the model. 
Notice that both the bispectrum and the trispectrum can be written as the sum of 
a purely isotropic and an anisotropic parts. The orders of magnitude of these two parts 
can, for instance, be read from Table 2 for the trispectrum: each one among t^l'^, 
and Tj^l provide the order of magnitude of the level of both their isotropic and 
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Figure 12. Plot of f{9,S) = [{Bc_{e,5,x* ,X2,x-i)xlxlkl)/{glH'^n?N\)] evaluated at 
{x* = 1,X2— 0.9, X3 = 0.1) in a sample angular configuration (see Eqs. (11.3) through 
(11.13)). 




Figure 13. Plot of the anisotropic part of the trispectrum from the contribution 
due to vector-exchange diagrams in a sample angular configuration (see Eqs. (11-14) 
through (11.17)). 

anisotropic contributions, which are therefore comparable; t^l^ instead quantifies a 
purely anisotropic contribution which, as discussed in Sec. 10, can be comparable to 
the other three parts, if not the dominant one. A similar discussion applies to the 
bispectrum (see f^j^ and f^^ in Table 1). We can then conclude that, for the three and 
for the four point function, there is room in the parameter space of the theory for the 
anisotropic contributions to be as large as, or even larger than, the isotropic ones. 
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12. Trispectrum for /(r) models of gauge interactions 



We will now show that it is quite straightforward to extend the calculations we performed 
for / = 1 to cases where / is not a constant. One interesting model is the one 
studied in [149] and also recently discussed in [90] , where the field is effectively massless 
(mo = ^ = 0) so the action (8.11) for the gauge field becomes 



S = J d'x^g ^-lf\<p)g^'^gf"^F;,F;, + 



:i2.11 



where again F^^ = d^B^ - d^B'^ + g^e^^^BlBl. 

Let us introduce the fields A'^ and A", related by the equations = fB°- = aA^. The 
A° are the physical fields. 

We can expand the perturbations of A" in terms of creation and annihilation operators 
in the usual way 

SAnv,x)= [j% i: [emapAl{rj,q) + h.c.]. (12.2) 

If / = foa", with a equal either to 1 or —2 (/o is a constant), it is possible to prove 
[90] that the equation of motion for 6A1 is the same as the one for SB^, where by SB^ 
we mean the transverse mode function in Eq. (8.14). This is equivalent to saying that, 
under the assumption a = 1, —2, the physical gauge fields are governed by the same 
equation of motion as a light scalar field in a de Sitter space and so they generate a scale 
invariant power spectrum. Let us sketch the calculation of the trispectrum in this theory. 

The general expressions of the Abelian terms still hold, except that the power spectrum 
in Eq. (9.1) reduces to 

Pt' = T:^P+. (12.3) 
having gauged the longitudinal modes away. 

Let us now have a look at the non- Abelian part. First of all, we need to set n{x*) = 0. 
The interaction Hamiltonian to third and fourth order are the same as in Eqs. (9.38) 
and (9.39), but with extra /^(r) factors. The anisotropy coefficients such as /„ and L„ 
that survive after setting the longitudinal mode to zero do not change. On the other 
hand, the wavefunctions for the gauge fields are now given by SB = 6A/f = SB'^/foa"'. 
The new trispectrum therefore differs because of extra scale factors inside and outside 
the time integrals, which in general imply a different power of if* in the final results 
and a different momentum dependence in the isotropic part of the expressions. 
Let us take into account the trispectrum. The three non- Abelian vector contributions 
from Eq. (9.34) in the / = 1 case can be schematically written as follows 



{6B^)une3 - [SB'} J dr] {6By = [6B' } J dr] [6B' } , (12.4) 
^ (SB)^ J dr^' {6Bf J dr^" {6B)' 

= {SB^yj dv' {SB^^y I dr^" (55^)', (12.5) 
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Figure 14. Plots of the isotropic functions of some of the vector-exchange 
contributions in the equilateral configuration, for the / ~ model. In this and 
in the next figures, "v. e. (I, II, III)" represent the isotropic functions associated with the 
very last term in square brackets in Eq. (12.13); "v.e.(new)" represents the isotropic 
function associated with the A:ii44, fc2244, kn^T, and ^2233 terms in the second line of 
Eq. (12.13). 

{SB'),,,, ^ i6B)'J dr,{5Bf=(5B^)' j dr^{5B^)\ (12.6) 

where we indicate with the subscript line?) the contribution from the third 
hne of Eq. (9.34) and with v.e. and c.e. respectively the vector exchange and 
the contact interaction contributions to the first hne of (9.34). (^-B^) = 
{6 Bf {ki)6 B^ {k2)S Bl{k3)6 Bf (ki)) and we have omitted all the gauge and vector indices, 
as well as complex conjugate symbols, for simplicity. Let us see now how Eqs. (12.4)- 
(12.6) change if / = /qo" (a = 1, -2) 

{6B%^,, ^ {6Bf I dvf {6Bf ^ ( V / —r)^ (12.7) 




{SB%,, ^ {6BY / rfr/ r {6BY / rfr/ f {6B) 




:i2.8) 



{5B\,,, ^ {6Bfj dvf {SBf ^ (^^j / dr^a^- (^^) . (12.9) 
Using a = {—Hr])^^ in the previous equations, we get 
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Figure 15. Plots of the isotropic functions of some of the vector-exchange 
contributions in the speciahzed planar configuration (plus sign), for the / ~ model. 



{5B%^,, ^ Ht'^ (5B^{-r^*Y)' j dv (55^)' (-r/)", (12.10) 

X I drj" [SB^y i-rj"r, (12.11) 

{SB'),.,. ^ i/f {6B^i-vTyj dv {SB^y i'V?"- (12.12) 

Let us now consider more in details the a = —2 case for contact-interaction and vector- 
exchange contributions. The expressions for the anisotropy coefficients are respectively 
given by 

rpEEEEE ^ kiksih -k^-ki- kf^h ■ h2)[5ijSki - Sijkkikii - Sijhakia - Skiki2kj2 - hikikji 
+ Sijhskuh ■ h + 5kikiikj2ki ■ h + kim + A;2244 + Ajuss + A;2233 - k22Mh ■ h 
- kiiMh ■ h - ki2iiki ■ h + ki233ki ■ h + ^1234^1 ■ hh ■ h] (12.13) 

and by Eq. (9.73), for one of the possible permutations. These expressions are more 
complicated w.r.t T/j"/ in Eq. (9.63) and T^^op in Eq. (9.88) for the longitudinal modes. 
As a result, when studying the shape of the trispectrum, for the isotropic functions 
appearing in it, several diagrams need to be taken into account, one for each term in 
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Figure 16. Plots of the isotropic functions of some of the vector-exchange 
contributions in the speciahzed planar configuration (minus sign), for the / ~ 
model. 

T^ki^^ and T^^i^^^. For comparison with the / = 1 case, we plotted the isotropic 
functions associated with the very last term in square brackets in Eq. (12.13) (see 
"— w.e(/)", "— w.e(//)" and "— w.e(JJ/)" in Figs. 14, 15 and 16). By comparing these 
plots with the ones in Figs. 8, 9 and 10, it is evident that they have very similar shapes. 
On the other hand, when we consider the isotropic functions associated with terms that 
are not present in the / = 1 case, several differences arise in the plots; we provide 
a sample in Figs. 14, 15 and 16 with the "w.e.(?2ety)" plots, which represent isotropic 
functions associated with the A;ii44, A;2244; ^1133 and A;2233 terms in the second line of 
Eq. (12.13). We verified that similar observations can be made concerning the shapes 
of the contact-interaction contributions. 
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13. Conclusions 



Cosmology has entered what it is called its "precision era": over the past few years, 
observations of the CMB and of other cosmological probes have been performed which 
have greatly improved the previous bounds on some of the fundamental parameters 
characterizing the early Universe cosmology. These bounds are expected to become a 
lot stricter with the advent of new experiments, such as the ongoing Planck satellite 
mission. The theoretical search for models of the early Universe has consequently been 
concentrated on, among other things, trying to produce more and more accurate pre- 
dictions. 

In this optics, alternative models to the basic single-field slow-roll inflationary scenario 
have been and are being investigated, with a strong focus on the computation of higher 
order (three and four-point) correlation functions and on higher order (beyond tree- level) 
corrections. All of these "higher order" predictions are expected to be particularly re- 
vealing of the early Universe physics, given that they are an indication of the type of 
fields and interactions populating the specific inflationary scenarios. 

In this thesis I have presented an overview of the existing results concerning some of 
these "higher order" predictions; in particular, I focused on our works on loop corrections 
in inflationary scalar field models and on the primordial non-Gaussianity and statistical 
anisotropy predictions in some inflationary models populated with non-Abelian vector 
fields. 

Loop corrections to the power spectrum of the comoving curvature fluctuation ( in 
single-field inflation arise both from the inflaton self-interactions and from the coupling 
of the scalar field with gravitons. We have calculated the corrections from tensor loops, 
previously neglected in the literature for simplicity reasons. It turnes out that one- 
loop corrections from tensor-scalar interactions are of the same order of magnitude as 
those arising from scalar self-interactions, therefore they cannot be neglected in a self- 
consistent calculation. 

One loop corrections have been found to be suppressed by an {H/mpY factor compared 
to the tree- level result; they exhibit a slightly different dependence from the external 
momentum, because of the presence of a logarithmic factor (which does not spoil scale 
invariance) and they turn out vary with time at most as fast as the logarithm of the 
scale factor. One-loop diagrams generally present both infrared logarithmic and ultravi- 
olet power-law (in addition to logarithmic) divergences. The ultraviolet divergences can 
be treated using ordinary regularization and renormalization techniques as in flat-space 
quantum field theory. The infrared divergences are "cured" by introducing an infrared 
cutoff represented by the smallest observable physical mode, i.e. by considering a finite 
space, a sort of "box" of observation, with a size equal to the current horizon length. 
As expected, tensor modes provide, as well as the scalar modes, observably small cor- 
rections at one loop level, as long as the basic inflationary scenario is concerned. Their 
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computation is nevertheless instructive and provides a path that can be followed in the 
investigation of "non-standard" models of inflation. More appealing results are for in- 
stance expected in inflationary scenarios described by non-canonical Lagrangians, some 
of which, P{X, 0) models, have been described and analysed in this thesis. 

Motivated by an interest in models that combine non-Gaussianity and statistical 
anisotropy predictions for the CMB fluctuations, we have worked on theories of in- 
flation where primordial vector fields effectively participate in the production of the ( 
perturbations, eventually focusing on some SU{2) vector fields models. The two, three 
and four point correlation functions in these models result as the sum of scalar and vec- 
tor contributions. The latter are of two kinds, "Abelian" (i.e. arising from the zeroth 
order terms in the Schwinger-Keldysh expansion) and "non-Abelian" (i.e. originating 
from the self-interactions of the vector fields). The bispectrum and the trispectrum final 
results are presented as a sum of products of isotropic functions of the momenta (F„ 
and Gn in the text) multiplied by anisotropy coefficients (which we indicated by and 
Ln) that depend on the angles among all gauge and wave vectors. 
The amplitude of non-Gaussianity has been evaluated through the parameters /tvl and 
tml] in particular we have discussed the dependence of these functions from the non- 
angular parameters of the theory. We have provided the comparisons among the different 
(scalar versus vector, Abelian versus non-Abelian) contributions to /nl and tnl, notic- 
ing that any one of them can be the dominant contribution depending on the selected 
region of parameter space. In particular, we have stressed how the anisotropic contri- 
butions to the bispectrum and the trispectrum can overcome the isotropic parts. An 
interesting feature of these models is that the bispectrum and the trispectrum depend 
on the same set of parameters and their amplitudes are therefore strictly related to one 
another. 

We have presented the shapes of both the bispectrum and the trispectrum. The isotropic 
functions appearing in their final expressions have been analyzed separately from their 
anisotropy coefficients. The bispectrum isotropic functions have been found to prefer- 
ably show a local shape. The trispectrum ones have been plotted selecting equilateral 
and specialized planar configurations. Finally, the full expressions (i.e. completed by 
the anisotropy coefficients) of bispectrum and trispectrum have been presented in spe- 
cific momenta configuration, in order to provide a hint of the modulation of shapes and 
amplitudes operated by anisotropy. 

In our view, the most promising features of these models consists in the possibility 
of providing non-Gaussianity and statistical anisotropy predictions that are related to 
one another because of the fact that they share the same underlying theory. Models 
that combine both types of predictions could be more easily testable and, from non- 
Gaussianity measurement, more stringent statistical anisotropy predictions could be 
produced or viceversa. 

Inflationary models that do not spoil the current agreement with experimental data 
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constitute a huge variety but many of them have very distinctive features that might 
be confirmed or ruled out by observations. Both the the non-Gaussianity and statisti- 
cal anisotropy predictions and the nature and the amount of higher corrections to the 
cosmological correlation functions can be ranked among these distinctive features and 
certainly deserve further investigations, looking forward for a confront with new and 
promising experimental data. 

While completing this thesis, the Wilkinson Microwave Anisotropy Probe (WMAP) 
team published its seven-year data analyses [153, 154]. The new bounds on non- 
Gaussianity from this study are given by -10 < /X < 74 and -214 < fj^l < 266, at 
95% CL. 
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15. Appendices 



15.1. Computation of the interaction Hamiltonian 

The propagator of two fields 0i and 02 is defined by (see for example [151]) 

(0102) = / /^0£'ne'/'^*"("'^-^, (15.1.1) 

where 11 is the momentum conjugate to and H is the Hamiltonian density. If H is 
quadratic in 11, as it happens for example in fiat space-time for a field governed by a 
Lagrangian L = Jd^x (|9^9^0 - 1^(0)) the square in the exponent can be completed 
and the integral in 11 evaluated and all is left is 

(0102) = J D0e^^ (15.1.2) 

So, if an interaction term with time derivatives appears in the Lagrangian, 11 and are 
independent fields in the path integral. This will provide some extra vertices that need 
to be accounted for in the Feynman diagrams. We sketch a derivation of these extra 
vertices. It will turn out to be similar to what is done in [108], although complicated 
by the presence of gravitons. To keep the calculations easier we will at first ignore all 
spatial derivative and tensor indices, this will also make the notation simpler. Also, we 
will momentarily ignore all numerical real coefficients; it is instead very important to 
keep track of imaginary coefficients, time derivatives and powers of the scale factor a, 
so we will make sure they are all accounted for in our analysis. 
The total action is 5* = J dr] (L^ + L^), where 

= a^7'^ + r^7'50'^ + r^7'^ + r^6(j)''y''^ + T^^dcj)^' + T^^^'"^ 
+ A0^^50 + A(^(^^7 + A^^^7 , 

= a^S(j)'^ + Ti6(j)'^ + 12(50'^ + io6(f)'^ + A50', (15.1.3) 
where /' = df /drj and where we define 



~ 07, 


(15.1.4 


~ a50. 


(15.1.5 




(15.1.6 


F2 ~ a^5(j)^ , 


(15.1.7 


V^^ ~ a^(507. 


(15.1.8 


i ~ a 7 , 


(15.1.9 


UJ ~ a(50. 


(15.1.10 


A</,^</, ^ aS(f)^, 


(15.1.11 


Xcjyyy ~ a(507^. 


(15.1.12 


A00T- ~ 07^0^, 


(15.1.13 




(15.1.14 
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Notice that in equations (15.1.4) through (15.1.14) we use the equivalence symbol 
meaning that we skeep details about integrations in momenta and real coefficients. 
The conjugate momenta are 



07 



07+ A^-^-y + A(^(^^ + r. 



f + r^7'50' + r^^7'' + r^7'', (15.1.15) 



5L 



5 50 



r^7 



r-y</,7 + A, 



+ A000 + u5(j)llh.l.l<o) 



We solve perturbatively the equations (15.1.15) and (15.1.16) in order to derive 7' and 
5(h to fourth order 



7 



-i4> 



;i5.1.17) 



5<p 



■ 7 

-4r 



(15.1.18) 
where 



The next steps are: derive the hamiltonian H = IP^y' + Il'^6(f) — L (^7, 6(f), 7 
we need to plug in the solution (15.1.17) and (15.1.18) for 7' and 6(f) ; construct the 
action as S" = 6*0 + Su, where Sq = J dr] {L^ + -^0) and Su includes the terms that 
depend on the conjugate momenta of the fields (a change of variables similar to the one 
that Seery performs in [108] over the conjugate momenta will also be necessary). 
Let's consider the vertices in Su that are involved in the corrections to the one loop 
point function for the scalar field 

Su^ J d7]' [a-^Ti6(f)'U^ + a-^riH'^n'^ + a-'^Ta'^l^ + a-^FsH'^n'^ 

'n^n^ + a-^u6(f)'U'^U't' 



+ 0-^^50 n^n^ + a-^w50 mwp^ . (15.1.19) 

The first three vertices belong to the third order part of the action; a^*Ti6(p U'^ and 
a~^ri7'n'^, provide a correction to the two point function at one loop with two vertices. 
Because of the presence of Fi which involves a factor of 0, it is subleading in slow roll 
order w.r.t. the corrections coming from fourth order vertices. We will therefore neglect 
these diagrams. The same applies to the second vertex, a~^rin'^n''^, although this may 
contribute to correcting the one point function 

(50,Kr/*)) D Ci-^v^ J d\h{q) f^dr]'6{0) [l - tk^') e''"' +c.c.,{15.1.20) 

where Ci is a numerical real coefficient and 5(0) is the Dirac delta function deriving 
from the propagator of the II'^^'s and /i is a scalar function of the internal momentum. 
The main contribution to the integral is due to times around horizon crossing since at 
early times the rotation to imaginary plane of the contour integral makes the exponent 



69 



decrease rapidly to zero and moreover rj* was chosen to be just a few e-folding after 
horizon crossing. Also, since the integrand function goes to zero as rj' approaches zero, 
we get a good approximation of this integral taking the upper limit 77* —t- 0. The result 
is purely imaginary and it cancels out with its complex conjugate. 
Let us now move to the fourth order vertices. From a~^r2n'^n'''' we have 

(%-;(r?*)50,-^^(r/*)) D J d\Uq) ^^r/'5(0)(l - ikri)\^^^'^' + c.c. (15.1.21) 

The same consideration as in (15.1.20) apply to the integral above, which gives a zero 
contribution, as well as the following diagrams (corresponding to the last two vertices 
in (15.1.19)) 

{HuSvlH^M)) ^ J d\h{q) dr^'6{0U' (1 - zW) e'^'"^' 

+ / d^qUiq) r dr]'S{0)ri'^ (l - ikr]') e^^^"' + c.c. (15.1.22) 



15.2. Study of leading slow roll order vertices in the fourth order action 

We are interested in computing the correlators just a few e-foldings after the scales we 
consider cross the horizon, so we can assume that the slow roll parameters remain small 
and can be treated as constants during this length of time. It is then correct to limit 
our interest to the leading order slow-roll contribution to the action . 



Let us start from the study of the slow-roll order of the fluctuations derived as solution 
to the constraint equations: 

where S'[50, 7] is a linear function of 5(f) and/or its derivatives and a linear function of 
7 and/or its derivatives, R[5(j)^] is a quadratic function of 50 and/or its derivatives and 
so on. Notice that the first order fluctuations are subleading (~ ^/e) w.r.t. the second 
order ones (~ e"). This criterium allows a suppression of a large number of terms in 
the 4th order action based on keeping the leading order (i.e. ~ 5°) terms only 



I dtd^x 



- y-^liklkjdjS(j)diS(j) + a2diS(f)di' 

+ \ (l^kiabdbjak - \iabl3kdklab 



dk02iabdb1ak + -;pabdk1abdk02 
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It can be easily shown that the terms in the action that do not contain the gravitons 
reproduce the ones in equation (37) of [152]. The contribution to the power spectrum 
due to these vertices has been calculated by these authors, but only for the scalar part. 
We will then focus on all the tensor contributions from these and from the remaining 
terms. Interaction vertices with both two and four tensor fluctuations will be obtained 
once the expressions for 02, ^2 and Pj are plugged in the action. The terms in the action 
that we need for constructing Feynman diagrams with one loop of gravitons are 



d xdt 



J^2 



4a2 



dj6(j)di6(j)'yik'ykj 



:i5.2.11 



2a2 

Let's plug the expressions for [ij and 6*2 into (15.2.1) considering the terms with two 
gravitons. The result is an ensemble of vertices which can in principle contribute 
to the one loop corrections to the power spectrum of the scalar field. Apart from 
djS(l)diS(j)'yik'~fkj, all of the other terms contain time derivatives of one, two or three of 
the four fields 

Pjd'^Pj D a'^[d-^(^d,ndjS^dmS(f) - d'^6^dj6(f) + dj6^d'^6(j) 
- dmHdmdjdcp) {iikdakj - liidiikj) 



(15.2.2) 
(15.2.3) 

(15.2.4) 

)h5.2.5) 



+ d~^(d,M(f)dm6(l) - d'^5(l)dM + dMd'^5(t) 



dmS(j)d^dj5(j))d iiikdakj - ludiikj) 



S(j)djS(j)dj92 D 



16H ^ 



2^2 ^iTig^aTii} ~l~ llji'lj 



D 



iikdi'jkj - likdiikj 



iabdalbkdk02 D iabdalbkTTrd ^d, 



AH 



;diS(j)dii 



iabdalbkl^k 3 iabdalbk'^a^d' 



dmdk 



+ dk5(l)d^ 



draHdmdk 



iabdklabPk 3 iabdklabd ^ \dmdk5(t)d„ 
- draS4>dmdkS 



iabdklabdk02 D iabdklab 



- d^S(j)dkS(j) + dkS(l)d'^ 



+ dj5(j)dj5 



■2 1 

b 7TdiS(j)di: 



(15.2.6) 



(15.2.7) 



(15.2.8) 



We will now prove that the vertices that include time derivatives do not actually 
contribute to the two point function. First of all notice that the tensor fields carry 
polarization tensors e^j with the property q^eij = and are always contracted with 
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other tensor fields in the calculations; this implies that, if a partial derivative index is 
contracted with a tensor index, that diagram will be zero. Based on this observation, 
we can ignore several of the vertices with time derivatives. We are eventually left with 
only two of them, that we will call Vi, V2 and V3 

~ d, (6^dj6<p) (dalbcda^bc) , (15.2.9) 
V2 ~ d, {6<j)dj64>) [iabiah] , (15.2.10) 

~ ial.dk-1ah{Pk + 0^62). (15.2.11) 

where (15.2.11) is given by the sum of (15.2.7) and (15.2.8). Notice that the fields 
need to be contracted between each other and that I] a, a' ^ig ^ig =constant [120]; the 
derivatives of 5(f) contract with derivatives of 7, so this produces k ■ q factors. Therefore 
we have 

{5ct),^S'^*)5ct)^^{r^*))y,^y,^i5^'\h + k,)Ht [ ^f,{q^)k.qr rfV/2(V) + c.c.(15.2.12) 

J q J —00 

where fi{q^) and f2{f]) ^i-re some functions of q^ and 77'. This contribution is evidently 
zero for symmetry reasons. 



15.3. Complete expressions of one-loop two-vertex diagrams to leading order 

In the following we provide the explicit expression for Eqs. (5.26) and (5.30). Eq. (5.26) 

reads as 



{H^S'^n^^kM))^,^.) = + k2)^{a, Hk) + Hk£) + as 



where 



k^ 



ai= - —(5 + 5x*^ + 2x*^). 
15 ^ ^ ■ 



02 



15x 



*2 



2 + ( 5 — 8(TcO"c + 47r(Ts — Sa^as ]x*'^ — Siirac + 2a(.(J, 



2ac(Tc)x*^ + (1 + ^cic<7c — 47r(Ts + Sas(Ts)x 



*4 



(15.3.1) 
(15.3.2) 

;i5.3.3) 



as 



1800x*2 
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3120 + 15136ac'5-c - 4507rVc(5-c - 75687r(T, 



2257rV, + 15136a,a,)x*2 - (l51367r(Tc - 4507rVc + 30272^^^, 



gOOvr^S-ea, - 30272(T,a, + QOOvrVca, ]x 



„*3 



672 - 15136(7,5-, 



+ 4507rVcae + 75687ra, - 2257rV, - 15136(7,5-, + 4507rV,5jx 



*4 



- 208x* 



(15.3.4) 



and p is a constant left over from renormalization of ultraviolet divergences. We have 
defined 

(7o = sin2x*, 
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= cos2x*, 
as = Si(2x*), 
a, = Ci(2x*), 

where Si and Ci stand for the sine-integral and the cosine-integral functions, i.e. 

Si(.) = / '^dt, 
cos(t) - 1 



Ci(x) = r -dt + ln(x) + 7, 

Jo t 



with 7 indicating the Euler Gamma function. 
The expression for Eq. (5.30) is 

(50,-;(r/*)50,-^^(r/*))f,^,2.) = ^5^'\ki + h)^ [c, + c,Hki)) , (15.3.5) 

where 

4/4 \ 

C2 = T7 -l^ + l • (15.3.7) 



15 Vx*^ 

Finally the quantity /a appearing in Eq. (6.1) is given by 

/3 = 03 + Ci + a', (15.3.8 
where a' = 2 (1 + x*^/3) a from Eq. (5.13). 



15.4- Background and first order perturbation equations for the gauge fields 



The equations of motion for the gauge fields and for a Lagrangian as in Eq. (8.11) with 
/ = 1 have been completely derived for the U{1) case in [135]. We are going to carry 
out a similar calculation for the SU{2) case 
1 



-A, 



-9 



a/3 



+g,e-''^g^^g'^F\f^'Bl + gy'^e'''^'g'''^g'PBtBlB'; = 

where F^;^^)" = d^B^ - d^B^. 

The 1^ = component of the equations of motion is 

djB] - djdjB^ + a^M^B^ + g^e"'''' 



(15.4.1) 



2BpjB'Q - B]B''. 



(15.4.2) 



where B^ = B'^{x,t). 

Let us now move to the spatial {u = i) part of (15.4.1) 

^d.d.B'i + M^BI - d,B^ - Hd,B^ + ^d.djB'; 

r-abc 



Bt + HB' 



+9ce 



abc 



9. 



C r, 



(a,S))5f + Bp,Bl 
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_abc 



(diBQ^BQ - B-Bq 



, 2 abc bb'c' 



Tyb' rye' 



r-.abc 
He ^abc^bb'c' 



B^B^B^ 



- d,B']B[ 



(15.4.3) 



If we contract Eq.(15.4.1) with d,^, we get the integrabihty condition 

{aMfB^ - M^diB^ + 3H[dAB^ - diBf) + g,e'''''[2H [diB^B^, + B'^diB', + B!>B]) - djB'^B'j 



-B^B^ + d.B'oB- - d^B'^B"^ - O.BfB, 
-didjB'.B;;)] + gye''''' [a^B'^B'^ B, 



B^d^B^ + diBpjB^ + BpjdiBI + d''B]B^ 



r>b r>c 
BqBq Bq 



ryb rye 
-Bi Bq 



BtBl Bt 



Bi Bq Bi 



diB^B^ Bq 



I RC R" RC 
c b' c 

B^diBi Bq 



f 2HB^,B', Bt 
B^Bid.Bi 



+— diB'B' B' + B^B' B^ + B'iB d^B 



(15.4.4) 



which reduces to Eq.(7) of [135] in the Abehan case. 
Combining Eq. (15.4.4) with Eq. (15.4.2) we get 

{aMfBS - M^diBf + 3H[a^M^B^ + - {djBfjB'^ - 2B^jdjB'o - B'^B^ + g.e'''''' B^^B^^ B^' 



+9ce 



abc 



2H(diB'lBl + B%Bl + Bpf) - djB'^B 
~Bp^ + d,BtB'^ - d^B^B', - d.B^B-, + \ 



d,d,B\B^A\ + gy^^e'' ^ [a^B^B^ B^ + B^B^ B^ + B^B^ B^ ) + 2HB^^Bl Bl 



^2^abchb 



c ub 



RC Rb RC RC RC DC ryb rye o jyc ryb rye Rcn jyb rye JDC ryb o jy 
-Dj-Do-Dj --Dj-Do-Di -^i^O^i -Oi^o^i^O - ^oOiBi Bq -Bf^BiOiB^ 



yc ryb 



c ryb 



yc ryb 



+^(diB';B' B' + BmB' B'^ + B'^B' d^B 



c ryb 







Plugging this into Eq. (15.4.3) we get 



(15.4.5) 



5^ + HB'^^ - ^djdjB^ + M'Bl + 2HdnBl 



aM 



+g,e'''' 2H[d,BlBl + B^diBl + E^^^j - d^B^B: 



dr. 



2^n 



-3H[g, 



.abc 



djB'ABl - 2Bp,Bl - B'^B^ + g.e^'''''' B]B^^ B"^ 



yb IDC 



yb rye 



-B)B'^ + djBlB'^ - d^B^^B^ - d.B^B'^ 



B^d^Bt + diB^djBt + Bp.diB'r + a^Ej^l 



-dA,B\B^^ 



2 abc bb c 
iJc 



RC rO rc RC Rb' RC RC rO rc a RC rO rc rcq rO rc rc rO q rc 

-BiB^Bi -B^BqBi -B^B^Bi - OiB^B^ Bq - B^OiB^ Bq - B^Bi OiB^ 
+\{d,BpfB^' + B'^d^B^Bj + Bpfd,Bf^ 



a^lB^B^ B'q + B'qB^q B'q + B'qBI B, 



yc jyb 



2HB'rBl Bl 

yc jyb 



+g,e^"^ HB'.B: + B',B: + B'.B^ 



+9c£ 



abc 



(dnB^^jB^ - B^^B^ 



-abc 



ic 9 



2abcWc' 

-g^e e 



RC r>V rye' 
BqB^ Bq 



9^ abc bb'c' 



-abc 

''^[{d,B])Bl + B]d,B^^ 
d,,B])B^^-(d,Bi)B 

c ryb' 



0. 



(15.4.6) 
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Let us consider the background part of the vector fields, i.e. 
Eq.(15.4.2) 



abc 



TDO TDO TDC 

BjBq Bj 



0. 



0. Then from 



(15.4.7) 



Before proceeding with the derivation of the equations of motion for the background 
and the field perturbations, it is necessary to make some comments about Eqs. (15.4.6) 
and (15.4.7). One approximation that we have been using in the computation of the 
cosmological correlators is allowing the fields to undergo slow-roll during inflation. 
One possible way to achieve this is by restricting the parameter space of the background 
gauge fields through the request that their temporal components should be much smaller 
than the spatial ones, Bq ^ |i?j^|/a(t), b = 1,2,3, and, in addition to that, assuming 
Bq ~ Bq, b,c = 1, 2, 3. With these assumptions, the temporal component can be 
factored out in Eq. (15.4.7), using the approximation B^ ^ HB^ (valid in a slow-roll 
regime). A solution to (15.4.7) is then given by Bq = 0. Adopting this solution and 
plugging it in Eq (15.4.6), it is easy to show that a slow-roll equation of motion for the 
physical fields 

2 



Af + ^HAf 



follows from (15.4.6) if we set 







2H^ with H <^ml and if 




{A^f H 




- {A^f C0S2 


- {A^f C0S2 






[A^f 




{A^fA 




- {A^f C0S2 ^23 


- (A1)^C0S2 di2 






{A^f 





(A1)'C0S2 

I A" I and cos 6 ah 



(15.4.8) 

(15.4.9) 
(15.4.10) 
(15.4.11) 



A" ■ A\ a and 



(A2)^-(A2)^cos2^23 

are satisfied. In the equations above, we defined A" 
h running over the gauge indices. The quantities appearing on the right-hand sides of 
Eqs. (15.4.9) through (15.4.11) can be either large or small w.r.t. one, depending on the 
specific background configuration, i.e. on the moduli of the gauge fields and the angles 

Suppose now the conditions described above are all met, then from Eq (15.4.6), in terms 
of the comoving fields, we have 

Bf + HB^ + M'^B'^ = 0. (15.4.12) 

Let us now derive the equations for the perturbations. Eq. (15.4.2) becomes 



djSB] - d'^5B^ + a^M^SB^ + g^e 



abc 



d,5Bp^, 



2Bp,5Bl 



5B]B^ - B^6B^ 



+9ce 



cb c 



5B'^B'q B] 



B]5Bl B^ 



B'jB'q 6B^ 



;i5.4.13) 



Eq.(15.4.1) for the field perturbations gives 

SBt + HSB^ - \djdj6B^ + M^6B^ + \didj6BJ - Hd,6B^ - Hdi6B^ 
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9c_ abc 



{d,5B])Bt + B]d,6Bl 



9 c nbc 



d,6B'^)B^^-{d,5B':)B^ 



9^ abc bb'c' 

9 



6Bp^ B^ 



B'^dBl B] 



Bp'; 6B'j 



-gc£ 



abc 



-9c£ 



_abc 



H {Bl5Bt + 5BlBt) + 5BlBt + B^dB^ + 5BlBt 

BUB' 



d,5Bl5Bl - 5B\Bl 



\ il'i pO-bc pbb c 



5BlBt Bl + Bl5Bt B^ + B^B^ SB^ 
Finally from Eq. (15.4.6) we get 

6Bf + H5B^ - \d^SB^ + M^SB^ + 2Hdi6B^ + (~ g^terms) = 0. 



(15.4.14) 



(15.4.15) 



When calculating n-point functions for the gauge bosons, the eigenf unctions we need 
are provided by free-field solutions, i.e. by solutions of Eq. (15.4. 15) with gc being set to 
zero. This is exactly the Abelian limit, in fact in this case Eq. (15.4. 15) corresponds to 
(18) of [135] and can be decomposed into a transverse and a longitudinal part 

2' 



+ Hdo + 



A;' 
a 



6B' 







9^ 
•^0 



1 + 



2k' 



F + {aMY 



Hdo + M'+ - 



5B^^ = 



(15.4.16) 
(15.4.17) 



where the time derivatives are intended w.r.t. cosmic time. 



15.5. Calculation of the number of e-foldings of single- (scalar) field driven inflation in 
the presence of a vector multiplet 

Let us consider the Lagrangian in Eq.(8.11) with / = 1 and = ml — 2H^. Let us 
assume that the SU{2) gauge multiplet undergoes slow-roll as well as the scalar field 
but the latter provides the dominant part of the energy density of the universe. This 
last hypothesis is necessary in order to produce isotropic inflation (i.e. in order for the 
anisotropy in the expansion that the vector fields introduce to be negligible w.r.t. the 
isotropic contribution from the scalar field). The expression of the number of e-foldings 
is 

N = Nscalar + N,ector = Nscalar + -Xr E ■ A\ (15.5.1) 

The previous expression can be easily derived from the equations of motion of the system 
neglecting terms that are proportional to the SU (2) coupling constant gc and assuming 
slow-roll conditions for both the scalar the gauge fields. 
The starting point is represented by Einstein equations 

2 SttG . . , , 

H = — - — [Pscalar + Pvector) ■ (15.5.2) 

where we split the energy density into a scalar and a vector contribution. In slow-roll 
approximation, p scalar ~ ^(0)- Let us calculate Pvector- The energy momentum tensor 
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for the gauge bosons 

5L 

where, as a remainder, L = -{l/4:)g^"^g''f^F^^F^p + {W^ /2)g^'''B1^Bl. So we get 



p>a p>a 2 2 tt tt2 

Too + B,+—B,B,- -B, B,+—B, B, 

I 9c abc rya ryb tjc , 9c ^abc^ab c ryb rye nb tjc 
+ Bi^0^i+-^^ ^ ^O^i^O^i 

2 

+ -^e^^'e"^'"' B^^B^Bf B^ (15.5.4) 

where sums are taken over all repeated indices. Let us write this in terms of the physical 
fields 

Tivector i i i / /ta i /\a Aa\ i abc f rj AO, i Ab AC i yc abc ab c /tb /\c /tb ac 

00 — 2 ^ y^i^i ^ ^O^Q) ^ 9c£ l^-n + J ^O^i + ^ ^O^i^O^i 

2 

+ ^£'^'"=£'^'''^'AM-A^'Af (15.5.5) 

If we neglect the non-Abelian contribution and we set Aq = 0, we are left with the 
Abelian result [123] 

The equation of motion for the background vector multiplet A"- can be derived from 
Eq.(15.4.12) 

Af + 3HA^ + mlA'} = 0. (15.5.7) 

which is equal to the equation of a light scalar field of mass mo, if mo <^ H. If the 
conditions for accelerated expansions are met, Eq.(15.5.7) reduces to 

3HAf + mlA'^ ~ 0. (15.5.8) 

We are now ready to derive Eq. (15.5.1). Let us start from the definition of N and keep 
in mind Eq. (15.5.2), where we are assuming the existence of a scalar fields in de-Sitter 
with a separable potential governed by the usual (background) equation 

4) + 3H^ + V' =0 (15.5.9) 

and slowly rolling down their potential. Then we have 

Jt* Jt* H Jt* 3H Jt' 3H 

^ ft V((j)) dt , ft (ml\ AfAfdA^dt' 

= 87rG -i^^d^ + SnCy^ 15.5.10 
Jt* 3H dd) ^ ^ JtA 2 3H dA°- ^ ' 



where A" = A'' ■ A''. So 



aJA^it')\^) 3HA]A] 



Hf) 3H(j) 
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1 fHt) V 
nip J<i>(t*) V 

nip J4,{t*) V 



1 



E 



AM? 



/ 1 \ _ /•A'-W 



E 



(15.5.11) 



after using the slow-roll conditions. Eq. (15.5.1) is thus recovered. 
In the final expression for the bispectrum then we can substitute 

dN / 1 ~ 



" dAf 



2m? 



A? 



(15.5.12) 



where the derivatives are as usual calculated at the initial time rj*. 
The upper limits in integrals such as the ones in Eq. (15.5. 11) depend on the chosen 
path in field space and so they also depend on the initial field configuration. It is 
important to notice though, as also stated in [148] , that if the final time is chosen to be 
approaching (or later than) the end of infiation, the fields are supposed to have reached 
their equilibrium values and so becomes independent of the field values at the final 
time t. Eq. (15.5. 12) is thus recovered. 



15.6. Complete expressions for the functions appearing in the bispectrum from quartic 
interactions 



J ^aa'b'ac'e 

J-EEE = £ £ 



The anisotropy coefficients /„ in Eq. (9.49) are listed below 
6 (iV'^' ■ N'') (iV^' ■ 1") 
+ [n^' ■ A^) [( - 2 ■ iV"') {ks ■ N^') -2{ki- iV"') {ki ■ N"' 

"(2 (iV"' ■ iV^') [k ■ N'') {k ■ A')) + (2 ^ 1) + (2 ^ 3)" 
{k ■ A'[2 [k ■ N"') {h ■ N'') {k ■ N'') +2[h- N^') [k ■ N"') {h ■ N'') 
h ■ N'^') [h ■ N'') {k ■ N'^') ki-h- {h ■ N"') {k ■ N'') [k ■ N^') k ■ h]) 
+ (2 ^ 1) + (3 ^ 2)]] (15.6.1) 
ku = e-'^'e-'^ [( {k, . N-') {ks ■ N'') {k ■ N"') {k ■ k) {k ■ A') 

- [k ■ N'' ) [k ■ N'') [k ■ N'') [k ■ k) {k ■ A') ) + (1 ^ 3) + (2 ^ 3)] (15.6.2) 



+ 



T ^aa'b' ^ac'e 

J-llE = £ £ 



+ 



(iV"' ■ A') ( (^1 ■ iV'^') (k2 ■ N'') + {k ■ iV"') [k ■ N'') )k ■ k 
(2 [k ■ A') {k ■ iV"') [k ■ N'') [k ■ N'^') ) + (1 ^ 2)" 

[ {{k ■ N"') {k ■ N"') + {k ■ N"') {k ■ N^')) {k ■ N'') [k ■ A') k ■ k) 

+ (1 ^ 3) + (2 ^ 3)]] (15.6.3) 
Ieei = e«"'^'e«^'^[4 (iV^' ■ A") (k^ ■ N^') {k ■ ^'') 

( {k ■ N"') {k ■ A') {{k ■ N'^') {k ■ N^') + (^3 • iV"') {k ■ iv^')) k ■ k) 



+ 
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+ (2 ^ 1) + (2 o 3) 

'(2 (k2 ■ A') [k ■ N"') {k ■ N"') {k2 ■ N'')) + (1 ^ 2) + (2 ^ 3) + (1 ^ 3) 
[n'' ■ h ■ k {{k ■ N"') {k ■ N^') + {k ■ iV^') (^1 ■ iV^'))) + (1 ^ 2) 



+ 



(15.6.4) 



where i ^ j means 'replace ki with kj\ whereas i ■v^ j means 'exchange ki with kj\ 



The isotropic functions F„ in (9.49) are given by 
1 



'FEE 



[Aeee + {Beee cos X* + Ceee sin x*) EiX* 



2Ak^klklklx*'^ 
Fill = n%x*)FEEE, 

FllE = n^{x*)FEEE, 

Feei = n^{x*)FEEE, 
where 

Aeee = kx*'^(^- k'^^kf + kl + k^ - ^kik2k^) - ^(^2/^3 + ^1^2 + ^1^:3) 

+ kik2k^{kl + kl + kl- A;2^3 - ^1^2 - kik3)x*'^^ 
Beee = {kf + kl + - A;^ + kik2k^x*^) 

Ceee = - k(kl + kl + ki)x*^{^ - k'^ + (^2^:3 + kik2 + kiks)x*^' 
In the previous equations we set k = ki + k2 + k^. 



(15.6.5) 

(15.6.6) 
(15.6.7) 
(15.6.8) 



(15.6.9) 
(15.6.10) 
(15.6.11) 



15.7. More details on computing the analytic expressions of vector- exchange diagrams 
We report the expressions of the functions A, B...P introduced in Eq. (9.57) 



{ki + k2) X* 



4k 

{ki + k2)x*' 
Ak 



A = i^-lGk'^ + kik2X*^j cos 

B = A[ki ^ ks, k2 ki] , 

C = Ak {ki + k2) X* cos 

D =C[ki^ ks, k2 ki] , 

E = (S/c^ {kf2 + ^3 + ki) — A;]^2^3^4^*^) cos 



ik {ki + k2) X* sin 



16k^ + kik2X*'^ ] sin 



{k^2 + h + ki)x* 



'{ki + k2)x* 
4k 

{ki + k2)x*' 
Ak 



Ak 



+ 2k {k^2 + ^3 + ^4)^ X* sin 
F = E[ks^ ki, ki k2] , 
G = 2k {k^2 + ki + k2f' x* cos 



{k^2 + h + ki) X* 
Ak 

'{k^2 + h + ^2)3;* 



Ak 



—8k {ki2 + ki + k2) — k^2^ik2X* ] sin 
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{kf2 + ki + ^2) X* 
Ak 



(15.7.1) 
(15.7.2) 
(15.7.3) 
(15.7.4) 

(15.7.5) 
(15.7.6) 

(15.7.7) 



M = (A;? + A;^2 + k%k2 + A;i2A;2 + ^2(^12 + ^2) + 

N = ki.ki k2] , 

P = L[ki ^ ks, k2 ki] . 



^2 



*2 • 

X St 



*2 • 



(A;i2 + ^1 + ^2) a;* 



Ak 

{k{2 + ki + /C2) 



4fc 



(15.7.8) 
(15.7.9) 



The anisotropy coefficients introduced in Eq. (9.66) have the following expressions 



ti = kiks (Jci ■ k 
t2 = kiki (Jci ■ k 
tg = k2k3 {k2 ■ k 
ti = k2ki (Jc2 ■ k 
= kik2 (ki ■ k 
tg = kik^ (^ki ■ k 
tj = k2k3 (k^, ■ k 
ts = k^ki (kj, ■ k 
tg = kik2 (ki ■ k 



(ki ■ k 

[k ■ k 

{h ■ k 
{ki ■ k 
{h ■ k 
{ki ■ k 
{ki ■ k 
ki ■ k 



[h ■ h 
[h ■ h 
[h ■ h 
[h ■ h 
{k2 ■ h 
(k2 ■ h 
[k2 ■ h 
k2 ■ k^ 
h ■ h 



^4 ' 



{h ■ ky. 

{ki ' ky 
h ■ k 
ki ■ k 

{k2 ■ ky 

(ki ■ k 

{k-k 



14 



(Jvi ■ k 

ho = kiks (k ■ ki) {k ■ k) [k ■ k) {k ■ k 
til = kihi (ki ■ ^^2) (^1 ■ ^2) [k ■ k) [k ■ ^12) ) 
ti2 = kiks (ki ■ ^12) (^1 ■ ^2) [k ■ k) [k ■ ^12) • 

Let us now list all the scalar products appearing in the equations above 



(15.7.10) 
(15.7.11) 
(15.7.12) 
(15.7.13) 
(15.7.14) 
(15.7.15) 
(15.7.16) 
(15.7.17) 
(15.7.18) 
(15.7.19) 
(15.7.20) 
(15.7.21) 



14 12 14 



14 2fcifc.-, 



ki ■ k 



12 



kl ' ky- 



13 



2kik^-,^ 



kl ■ k 



ko ■ k 



k'i +k'f-k'f 
1^1 1 



12 



2k2k 



k?. ■ k 



4^ 



13 



2^3 fc,-. 



ki ■ k 



14 



feZ. -k'i+k'i 
14 1 ' 4 

2^4^-1 4 



/C3 ■ k 



k '- ~\~ k '- 2Ajo 



12 



^2 ■ ki3 = ^-'^2k^3 ' 



k2- ki= 2fefcl ~ 



_ «f4+«i-2 «1 «3 ^^^4 
M Ai^2 — 2fc4 fci-;^ 



h. ■ h - — 14 12 1 3 4 



2 fc^: k -j"^^ 



ko. ■ k 



14 



3 14 2 
2fc2fc-|4 



kl ■ k 



2kik2 



kl ■ ks 



k A ~\~ kr\ k - k ^ 



~ I 1,2 _h.'2._h.'Z 

l. I. 12 4 j 

'^3 ^^4 - 2kiks 



2kik3 



kl ■ ki 



r i, . _ '^2 "3 "14 
1^3 l^li - 2k;k., 



2fcife4 



~ I L.2 _h.'2._h.-2. 

• ^^3 = ' 2fc2fc3 



where ki = \ki\, ki = ki/ki, k^j = ki + kj, kij = \ki + kj\, i and j running over the four 
external wave vectors. 



15.8. Complete expressions for functions appearing in point-interation diagrams 
We provide here the expressions for the coefficients appearing in Eq. (9.89) 
Qeeee =x*^[- k{ek\ + k^kl -kl + ekl + k^k^ki - 3k^klk3ki + k^klki 
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- kl + k^hkl + k^kl + k^k2{kl - kk^k^ - 3klk^ - ^k^kj + kl 

+ k\k^ + k^)) + kl{-l + k\k^ + k^)) - 3k'^kl{ksk4 + k2{k3 + k^)) 
+ kl{-l + k\k2 + A;3 + k^)) + k^ki{kl + A;| - ?>klki - ^.k^^kl + kl 

- 'iklik^ + A;4) + k'^{k2 + k3 + k^) - 3k2ikl + ^k^k^ + kl) - kik^ki 

+ k2{k3 + ki))))] + x*^[k'^{k\k2k3ki + k2kski{kl + ^3 - 3k2ksk4 + kl) 
+ klik^ki + k2{k:i + k^)) - 3kl{klkl + k2k3k^{ks + k^) + kl{kl + k^k^ 
+ kl)) + kiik^k^ki + k^2{h + k^) - ^klk^k^ik^ + k^) + fc3A;4(A;^ + 





~\~ k^^i^Iv^ ~\~ A^2^4 S/cg/c^ ~\~ ~\~ k"^k^(yk ~\~ ^4))))] 3x k-^k^k^k^^ 


{Ih.i 


U) 


^EEEE 


^ k ^k-^ ~\~ k'2 ^3 ~l~ *^ ; 


(15i 


^.2) 


Beeee 


= k^ -k'^ {k^ki + ^2 (A;3 + ki + ki (/ca + + ^4) a;*^ + kik2k3kiX*^y) , 


(15.^ 


U) 


Ceeee 


= /cx* (k^ - (/C2A;3/C4 + A;i {k^k^ + /ca (A;3 + k^))) x*^) , 


(15.^ 


U) 


Deeee 


= - kx* (k^ - {k2k3ki + kl {k^k^ + /ca (/i;3^4))) x*^) , 


(15i 


^.5) 


Eeeee 


= k^-P {k^ki + /cs (A;3 + A;4) + ki (/cs + A;3 + /^4)) a;*^ + kik2k3k4X*^, 


(15.^ 


^.6) 


Feeee 


^ /iy ~\~ k'2 ^3 ^4^ 


(15.^ 


^.7) 
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